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The  electrodynamic,  thermal,  and  energetic  character  of  stable  Sun-aligned  arcs  in  the  polar  cap  can 
be  meaningfully  diagnosed  by  an  incoherent  scatter  radar,  provided  a  suitable  observing  scheme  is 
selected.  We  report  here  such  measurements  of  two  intense  Sun-aligned  arcs.  The  two  arcs  were 
diagnosed  on  two  different  nights  (February  26  and  March  1.  1987)  using  the  Sondre  Stromfjord  radar 
as  a  stand-alone  diagnostic.  Repeatable  patterns  are  found  in  mesoscalc  area  (order  I03  km  by  103  km) 
maps  of  altitude  profiles  for  observed  electron  and  ion  gas  number  densities,  temperatures,  and 
lir.?-of-sigh'  velocities  "'r.J  (,iojc'-ted  mesoscale  area  maps  of  derived  electric  fields.  Pedersen  and  Hall 
conductivities  ( Ne .  Tt.,  T,.  V,  E.  Sw).  horizontal  and  field-aligned  currents,  joule  heating  rate,  and 
Poynting  flux.  We  confirm,  for  the  first  time  with  continuous  mesoscale  area  maps,  that  the  arcs  have 
the  anticipated  simple  arc  electrodynamics.  That  is,  the  visual  and  enhanced  ionization  signatures  of 
the  arc  are  produced  by  incoming  energetic  electrons  carrying  the  outgoing  current  from  the  electric 
field  convergence  in  the  arc.  Strong  electron  temperature  enhancements  (>2000  K)  are  found  as 
expected  within  the  sheets  of  ionizing  particle  precipitation.  Dawn  to  dusk  decreases  in  the 
antisunward  plasma  flow  of  order  1  kms'1,  across  order  100  km,  correspond  to  peak  electron  densities 
of  order  I05  cm-3  down  to  altitudes  as  low  as  120  km.  and  upward  currents  of  order  I  pA  m  These 
data  also  lead  to  important  implications  for  the  physics  of  polar  cap  arcs.  The  high-velocity 
(antisunward  flow  on  the  dawnside)  edge  of  the  arc  marks  the  location  of  strong  persistent  Joule 
heating  driven  by  downward  Poynting  flux.  There  is  a  channel  of  strongly  enhanced  ion  temperature 
(well  above  the  electron  temperature)  along  the  high-velocity  edge  of  the  arc.  quantitatively  accounted 
for  by  ion  frictional  heating  as  the  strong  electric  field  drags  the  ions  through  the  neutral  gas.  The 
deposition  rate  into  the  atmosphere  of  the  net  electromagnetic  energy  well  exceeds  the  net  particle 
energy  deposited  by  the  ionizing  energetic  electron  flux.  This  heating  is  a  substantial  source  of  heat 
into  the  polar  thermosphere.  As  estimated  by  heat  into  the  ions,  heat  lost  by  the  ions  to  the  neutrals, 
or  energy  available  from  the  Poynting  and  energetic  particle  flux,  several  ergs  cm  s' 1  are  deposited 
in  channels  of  order  100  km  for  good  fractions  of  an  hour  in  these  arcs:  this  contributes  to  resolving 
the  problem  of  the  missing  polar  thermosphere  heat  source.  Finally,  a  reasonably  simple  yet 
self-consistent,  accurate,  and  comprehensive  representation  of  stable  intense  Sun-aligned  arcs  is 
presented,  including  the  electrodynamic,  thermal,  and  energetic  character. 


1.  Introduction 

Morphological  studies  of  intense  polar  cap  aurora  began 
during  the  International  Geophysical  Year  (IGY).  Measure¬ 
ments  performed  during  this  year  at  Antartica  \  Weill.  1958; 
Denholm  and  Bond,  196 11  and  at  Resolute  Bay.  Canada 
[Davis,  19601  found  intense  polar  cap  arcs  oriented  in  the 
Sun-Earth  direction.  Later,  a  scanning  photometer  on  board 
ISIS  2  [Berkey  cl  at..  1976:  Ismail  et  at..  1977)  and  all-sky 
cameras  in  Greenland  [Lassen  and  Dani'dsen,  !9'7S]  re¬ 
vealed  that  the  occurrence  of  intense  polar  arcs  throughout 
the  polar  cap  correlates  well  with  northward  interplanetary 
magnetic  field  (IMF). 

New  initiatives  to  observe  the  intense  Sun-Earth  aligned 
arcs  involved  airborne  observatories  [Romick  and  Brown. 
1971 ;  Lather  and  Akasnfu,  1969|  and  also  particle  and  optical 
emission  detectors  on  Defense  Meteorological  Satellite  Pro¬ 
gram  <DMSP)  satellites  [Menu  and  Akasofu.  1976:  Akasofu. 
1976|.  The  new  instrumentation  served  to  indicate  that  the 
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polar  arcs  are  usually  produced  by  soft  (<1  keV)  electrons 
and  to  characterize  the  different  types,  locations  and  shapes 
of  intense  polar  cap  auroras. 

Statistical  and  case  studies  of  auroras  and  particle  precip¬ 
itation  in  the  polar  cap  indicated  that  intense  polar  cap 
auroras  preferentially  occur  in  the  morning  sector  I Giessen - 
hoven,  1982).  The  average  energy  of  intense  electron  precip¬ 
itation  in  the  polar  cap  was  found  to  be  between  50  and  300 
eV  [Hardy.  1984).  Burke  el  at.  (1982]  inferred  that  intense 
polar  cap  arcs  arc  the  optical  signature  of  upward  Birkciand 
currents  associated  with  corresponding  velocity  gradients. 

The  DE  1  spacecraft  provided  images  of  intense  Sun- 
aligned  arcs  extending  across  the  polar  cap.  The  transpolar 
auroras  together  with  the  auroral  oval  resembled  the  Greek 
letter  theta  | Frank  et  at..  1982.  1986J.  The  salient  features  of 
this  striking  auroral  pattern  are  the  presence  of  precipitating 
tons  accompanying  the  energetic  electrons  (>1  keV),  its 
location  in  a  region  of  sunward  convection,  and  the  possible 
existence  jf  four  cells  of  plasma  convection.  Burke  et  al. 
(1979!  have  also  presented  data  suggesting  four  convecting 
cells.  A  valuable  description  of  the  convection  for  different 
values  of  the  three  components  of  the  IMF  has  been  pre¬ 
sented  by  Reiff  and  Burch  1 19851-  Their  analysis  implies  a 
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four-cell  convection  pattern  for  strong  northward  IMF  and  a 
three-cell  convection  for  smaller  positive  values  of  B,. 

There  is  a  general  consensus  that  some  of  the  intense  polar 
cap  Sun-aligned  (S-A)  arcs  are  formed  in  a  region  of  closed 
field  lines  [ Akasofu  et  al..  1984:  Frank  et  al..  1986J  with 
precipitating  electrons  having  the  characteristics  of  plasma 
sheet  or  plasma  sheet  boundary  layer  electrons  [Peterson 
and  Shelley ,  1984;  Ohara  et  al.,  1988].  Kan  and  Burke  [1985] 
have  argued  that  the  less  intense  polar  arcs  possibly  occur  on 
open  field  lines.  Other  high-latitude  arcs  may  well  corre¬ 
spond  to  a  poleward  widened  auroral  oval  [ Meng ,  1981; 
Murphree  et  al.,  1982],  Data  presented  by  Huang  et  al. 
[1987]  collected  by  ISEE  1  during  excursions  in  the  magne¬ 
totail  suggest  filaments  of  plasma  sheet  plasma  extend  into 
the  lobes.  Even  though  the  filaments  have  been  observed  for 
any  orientation  of  Bz,  their  presence  reflects  the  discrete 
nature  of  the  plasma  source  and  suggests  a  causal  link 
between  the  IMF  and  the  formation  of  the  polar  cap  arcs. 

With  the  advent  of  intensified  imaging  photometers,  much 
more  common  weak  stable  Sun-aligned  arcs  were  found  in 
the  polar  cap  [Weber  and  Buchan ,  1981].  The  Sondrestrom 
incoherent  scatter  radar  was  also  used  to  study  the  morphol¬ 
ogy  of  these  weak  S-A  arcs.  Carlson  et  al.  [1984]  determined 
sheared  (antisunward  to  sunward  reversal)  flow  which  was 
coincident  with  bottomside  F  region  density  enhancements 
and  with  optical  emissions,  both  produced  by  the  precipitat¬ 
ing  particles.  The  plasma  velocity  gradients  coincided  with  a 
simple  arc  electrodynamics  interpretation.  Robinson  et  al. 
[1987]  obtained  the  field-aligned  currents  for  a  Sun-aligned 
arc  observed  near  local  midnight  and  merging  into  the 
auroral  oval.  Their  measurements  implied  an  upward  current 
on  the  western  edge  of  the  arc  and  a  downward  current  on 
the  eastern  side.  Mende  et  al.  [1988]  have  also  presented 
measurements  from  Sondrestrom.  Their  observations  of 
multiple  polar  cap  arcs  showed  a  preferential  antisunward 
flow  within  the  arcs. 

Carlson  et  al.  [1988]  have  presented  convection  patterns 
representing  northward  B-  conditions  which  were  obtained 
by  combining  the  optical  emissions  associated  with  S-A  arcs 
and  the  one-dimensional  velocity  measurements  of  DE  2. 
These  quantitatively  confirm  electron  impact  production  of 
the  arc,  and  local  simple  electrodynamics  for  the  arc.  Weber 
et  al.  [1989]  used  data  from  a  rocket  that  traversed  a  polar 
cap  arc  to  model  the  circuit  parameters  of  the  arc.  They  also 
showed  enhanced  electron  temperature  variation  across  the 
arc. 

The  purpose  of  this  paper  is  to  examine  more  comprehen¬ 
sive  diagnostics  and  analysis  of  the  electrodynamics,  the 
energy,  and  the  thermal  properties  of  a  few  Sun-aligned  arcs 
which  v.  ;r<*  •orated  in  the  polar  cap.  The  data  are  more 
comprehensive  by  virtue  of  time  continuous  measurement, 
over  a  wide  spatial  area,  of  a  significantly  more  complete  set 
of  geophysical  parameters.  The  arcs  to  be  described  in  this 
work  were  probed  with  the  Sondrestrom  radar  on  the 
evenings  of  February  26  and  March  1,  1987.  The  common 
signature  of  these  arcs  consists  in  that  they  were  produced 
by  electrons  with  energies  larger  than  I  kev,  as  deduced 
from  the  altitude  of  significant  E  region  ionization.  Although 
we  diagnosed  over  30  stable  Sun-aligned  arcs,  using  all-sky 
imaging  photometers  (ASIPs)  and/or  the  incoherent  scatter 
radar  (ISR)  during  this  campaign,  here  we  will  concentrate 
on  only  two  of  these  arcs.  The  arcs  we  concentrate  on  are 
the  two  most  intense,  for  optimum  signal  to  noise  ratio  or 


measurement  statistics  on  derived  parameters.  The  arc  of 
February  26  was  transpolar  as  observed  by  the  DMSP/F6 
satellite  [Niciejewski  et  al.,  1989].  This  stable  arc  had  an 
electron  density  above  103  cm*3  in  the  E  region.  On  March 
1,  1987  we  observed  three  polar  cap  arcs.  Two  of  them  were 
Sun-aligned  and  the  other  hook-shaped.  The  two  Sun- 
aligned  arcs  were  initially  separated  by  80  km  with  peak 
densities  located  at  substantially  diflerent  altitudes.  They 
subsequently  merged  into  a  single  broad  arc.  We  discuss  the 
properties  of  Sun-aligned  arcs  for  these  two  nights  to  dem¬ 
onstrate  repeatability  of  some  basic  characteristics. 

In  section  2  we  identify  a  data-taking  mode  designed  for 
such  studies,  and  present  Sun-aligned  arc  parameters  di¬ 
rectly  observable  by  this  technique.  In  section  3  we  eAtend 
analysis  of  the  ISR  data  to  present  important  derived  param¬ 
eters.  In  section  4  we  discuss  the  implications  of  these 
measurements,  develop  their  interpretation,  and  identify 
their  significance  to  understanding  of  the  electrodynamics, 
thermal  character,  and  energetics  of  the  arcs.  Finally,  in 
section  5  we  summarize  the  main  conclusions  from  this 
study. 

2.  Data-Observing  Mode  and  Direct 
Observables 

The  Sondre  Stromfjord  incoherent  scatter  radar 
(66.987°N.  50.949°W.  74°  Invariant  Latitude),  is  situated  100 
miles  (160  km)  inland  from  the  west  coast  of  Greenland.  Its 
location  at  a  very  high  magnetic  latitude  makes  this  station 
capable  of  probing  the  nighttime  polar  cap  and  examining  the 
dynamics  of  polar  Sun-aligned  arcs  [Cur/soti  et  al..  1984: 
Robinson  et  al..  1987;  Mende  et  al..  1988], 

While  it  is  difficult,  if  not  impossible,  to  recognize  Sun- 
aligned  arcs  with  conventional  ISR  observing  modes  in  the 
absence  of  coincident  optical  imaging  data,  it  is  possible  to 
design  a  radar  mode  that  can  unambiguously  recognize 
them.  We  describe  the  design  here,  illustrated  in  Figure  1. 
The  design  is  based  on  what  we  know  of  their  morphology 
from  optical  imagery.  The  essence  of  the  design  is  simple. 
One  wants  to  map  the  electron  density  profile  and  velocity 
field  over  an  area,  with  the  area  large  compared  to  the  arc. 
and  with  time  resolution  short  compared  to  the  time  for  the 
arc  to  pass  through  the  ISR  field  of  view.  Signal  to  noise 
considerations  determine  look  elevation  and  area  of  cover¬ 
age.  and  limit  azimuth  scan  rates  and  cycle  time  resolution. 

Variability  in  the  motion  and  lifetime  of  S-A  arcs  [  Weber 
et  al.,  1989]  also  imposes  some  restrictions.  For  this  reason 
we  implemented  the  azimuthal  sweeps  (AZ)  at  two  diflerent 
speeds.  One  AZ  scan  was  performed  at  a  rapid  scanning 
velocity,  with  the  purpose  of  obtaining  a  snapshot  of  the 
electron  density  inside  fast  moving  arcs.  Two  slow  AZ  scans 
'wept  90“  in  2  min.  probing  slower  moving  arcs  and  provid¬ 
ing  the  line  of  sight  (LOS)  .  .lociiy  and  the  number  density 
with  smaller  error  bars.  Good  horizontal  transarc  spatial 
resolution  was  also  achieved  (15  km  at  150  km  of  altitude 
when  pointed  across  the  arc.  and  ultimately  to  tE. 

I.  j-knt  raitai  beam  wiUth  smeared  over  the  integration  time 
when  pointed  along  the  arc).  These  measurements  were 
complemented  with  an  elevation  (EL)  scan  (perpendicular  to 
the  anticipated  arc  direction)  that  gives  the  dependence  of 
the  geophysical  parameters  on  altitude  and  distance  across 
the  arc.  The  scan  sequence  (shown  in  the  bottom  panel  of 
Figure  ll  was  slow  AZ-fast  AZ-slow  AZ-EL.  l  or  the  data 


Valladarfs  and  Carlson:  Polar  Arc  Electrodynamics/Energetics 


1381 


Fig.  I .  Schematic  representation  of  the  volume  scanned  by  the 
Sondrestrom  radar  during  the  EL  and  AZ  sweeps  of  February  26, 
1987.  The  bottom  panel  shows  the  sequence  followed  by  the 
antenna.  The  initial  azimuth  was  selected  in  order  to  align  node  3 
with  the  Sun-Earth  direction. 


presented  here,  one  complete  cycle  was  performed  in  8  min 
and  40  s.  Figures  2a-2h  present  the  data  corresponding  to 
two  consecutive  cycles  of  this  antenna  pattern. 

We  will  also  show  data  from  two  other  observing  modes. 
One.  tailored  to  a  specific  need,  looked  directly  up  the 
magnetic  field  line  so  as  to  give  true  instantaneous  altitude 
profiles  within  a  Sun-aligned  arc.  Such  time  continuous 
observations  during  the  full  passage  of  an  overhead  arc  thus 
allow  detailed  diagnosis  of  the  properties  of  a  transverse 
section  of  an  arc.  These  observations  are  presented  in  Plates 
la  and  lb. 

The  third  data  set  presented  here.  Figures  4a-4c,  for 
March  I,  1987.  again  shows  electron  number  density  and 
LOS  velocity.  The  antenna  pattern  consisted  of  a  full  360° 
AZ  scan  and  two  120°  EL  scans  embedded  at  planes  perpen¬ 
dicular  to  each  other.  By  chance,  a  Sun-aligned  arc  is  clearly 
apparent  150  to  200  km  duskward  (and  westward)  of  the  ISR 
on  the  Figure  4a  scan,  although  the  20-min  time  resolution  of 
this  stand:". I  Sondrestrom  observing  mode  is  too  slow  to 
track  or  in  general  unambiguously  recognize  an  arc.  The 
s0r«Hrpstrom  ISR  is  able  to  make  a  pioper  idemiueation  of 
polar  cap  arcs  based  upon  (I)  the  ability  of  the  radar  to 
measure  the  number  density  and  the  LOS  velocity  while 
doing  rapid  scans,  and  (2)  he  well-known  distribution  of 
density  gradients  and  velocity  shears  within  the  arcs  [Carl¬ 
son  et  ai..  I984|.  Moreover.  E  region  polar  cap  arcs  exhibit 
enhanced  densities  below  200  km.  Recall  that  chemical 
lifetimes  are  so  shoit  below  200  km  that  strong  A',,  enhance- 
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Fig.  2a.  Number  density  contours  and  line  of  sight  (LOS) 
velocity  measured  by  the  Sondrestrom  radar  on  February  26,  1987 
and  during  the  time  interval  from  0224  to  0226.  The  azimuth  was 
kept  fixed  at  23 1°  and  the  elevation  varied  from  25°  to  60°.  The  initial 
density  contour  corresponds  to  0.4  x  105  cm”5,  and  consecutive 
contours  increase  linearly  in  steps  of  0.4  x  10 5  cm-3.  The  LOS 
velocity  was  scaled  in  order  to  represent  a  velocity  equal  to  1  km 
s  _  1  as  a  segment  100  km  long.  The  dashed  line  vectors  indicate  LOS 
velocities  directed  toward  the  radar,  and  the  double  line  vectors 
characterize  plasma  motion  away  from  the  radar. 


Fig.  2 h.  Plasma  parameters  corresponding  to  the  AZ  scan  of 
0226-0228.  The  constant  elevation  was  35°  and  the  azimuth  changed 
from  23T  to  141°.  The  measured  velocity  and  density  have  been 
projected  to  a  horizontal  plane  using  the  same  scaling  factors 
employed  in  Figure  2a.  Geographic  north  is  at  the  top  and  west  is  on 
the  left  side  of  the  plot.  The  arrow  at  the  top  central  part  indicates 
that  during  the  experiment  the  direction  of  the  Sun  was  7°  west  of 
geographic  north. 
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Fig.  2c.  Plasma  density  for  February  2b.  1987  and  time  interval 
between  0228  and  0230.  The  antenna  scanned  180°  in  2  min  using  a 
fixed  elevation  of  35°.  The  scan  started  at  an  azimuth  of  141°  and 
ended  at  321°.  The  measured  LOS  velocity  is  not  shown  in  this  plot. 


ments  identify  local  production  of  ionization.  We  have  used 
the  density  and  LOS  velocity  mapping  capabilities  provided 
by  the  AZ  scans  to  properly  determine  the  occurrence  of 
arcs. 

The  densities  plotted  in  Figures  2-4  and  Plate  I  are  on  a 
linear  scale.  The  initial  value  and  the  interval  between 
consecutive  levels  are  both  equal  to  0.4  x  I05  cm  \  The 
arrows  superimposed  to  the  N,.  contours  represent  the  LOS 
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Fig  2d.  lop  view  of  the  density  and  the  I  .OS  velocity,  both 
projected  to  a  horizontal  plane.  The  scan  time  is  from  0230  up  to 
0212  This  plot  corresponds  to  the  northern  quadrant.  T  he  elevation 
is  and  the  scan  starts  at  421  of  azimuth  and  continues  until  231 
I  he  scaling  factor  is  the  same  as  m  previous  figures. 


Fig.  2c.  Similar  to  Figure  2u,  but  for  the  time  interval  0232-0234. 


velocities:  they  have  different  patterns  and  lengths  according 
to  the  sign  and  magnitude  of  the  LOS  velocity.  At  150  km  of 
altitude  the  density  and  velocity  statistical  uncertainties 
were  about  5C7  and  <100  m  s  respectively. 

2.  L  February  26,  1987  Scans 

Data  from  this  day  show  a  stable  Sun-aligned  polar  arc 
that  stayed  in  the  radar  field  of  view  for  40  min.  This  arc 
penetrated  to  the  E  layer  [N,.  was  enhanced  at  altitudes  as 
low  as  120  km),  and  was  detected  near  local  midnight.  Its 
southernmost  region  was  observed  to  be  merged  to  the 
poleward  edge  of  the  auroral  oval  during  the  early  part  of  the 
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F  ig  2 f.  Similar  to  Figure  2 h.  hut  con  expending  to  the  time 
between  0235  unit  0237. 
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Fig.  2 g.  Similar  to  Figure  2r.  but  for  time  0237. 


arc  transit.  At  the  time  of  the  experiment  the  background 
density  in  the  polar  cap  F  region  was  0.2  x  105  cm  ~3,  a  value 
well  below  the  lowest  density  contour  level  of  Figure  2.  The 
Pedersen  conductivity  peak  values  corresponding  to  each 
scan  are  presented  in  Table  1,  which  also  includes  the 
altitude,  distance  along  the  dawn-dusk  meridian,  and  dis¬ 
placement  velocity  of  both  maxima.  The  geomagnetic  Kp 
index  was  1+  and  the  solar  flux  (Sa)  73.9. 

The  data  of  Figure  2 a  were  obtained  between  0224  and 
0226  UT  when  the  S-A  arc,  as  defined  by  scans  of  Figures 
2b-2d.  was  located  to  the  west  of  the  radar  station.  The 
beam  scanned  through  most  of  the  arc,  and  also  probed  a 
duskward  region  beyond  the  arc  (left  side).  In  this  region 
beyond  the  duskward  and  equatorward  edge  of  the  arc  the 
LOS  velocity  was  directed  away  from  the  Sun  and  radar  (the 
small  electron  density  in  this  region  precludes  measurement 
of  a  precise  magnitude  of  the  velocity).  Within  the  arc.  at  F 
region  altitudes  (see  range  gates  2  and  3  placed  at  160  and  210 
km  of  altitude),  the  LOS  velocities  (looking  transverse  to  the 
arc)  were  less  than  100  ms"1.  That  the  density  contours 
trace  to  lower  altitudes  in  the  dawn  (right)  side  of  the  arc 
suggests  harder  precipitation  there. 

Figure  2b  presents  the  data  gathered  with  the  slow  AZ 
scan  of  0226.  Here,  the  Sun-aligned  arc  is  identified  by  the 
high  Ne  structure  seen  in  the  lower  south  to  southwest  part 
of  the  plot,  with  an  orientation  much  more  clearly  seen  in 
Figures  2d.  2/,  2g,  and  2 h.  The  southern  extreme  of  the  arc 
shows  (in  Figures  2b,  2c.  and  perhaps  2 f)  a  dawnside 
(eastward)  extension  with  an  electron  density  (about  >0.5  x 
I05  cm  ')  larger  than  the  background  A,,  at  altitudes  about 
200  km.  This  segment,  not  Sun-aligned,  but  more  auroral 
oval-aligned,  is  probably  associated  with  the  connection  of 
the  Sun-aligned  arc  to  the  poleward  edge  of  the  auroi.il  oval. 
The  northern  extreme  of  the  Sun-aligned  arc  A,,  contours 
mapped  here  in  Figure  2 h.  all  clearly  pointing  sunward,  is 
determined  by  the  simple  geometry  of  the  arc  intersect  with 
the  conical  radar  scan. 

Still  in  Figure  2b.  gates  2  and  3.  set  at  the  same  heights  of 


1 3X3 


the  previous  EL  scan,  indicate  a  general  antisunward  motion 
and  a  shear  inside  the  arc.  Outside  in  the  duskside  there  is  a 
small  sunward  flow;  conversely  the  dawnside  region  has  a 
large  antisunward  flow.  We  will  come  back  to  this  point  in 
section  3.  where  we  include  the  Earth  corotation  velocity 
and  compute  the  direction  of  the  plasma  flow  in  the  Earth- 
Sun  frame  of  reference. 

Figure  2c  shows  the  same  S-A  arc  during  the  fast  AZ.  180° 
scan  of  0228.  The  two  quasi-ellipsoids  that  are  observed 
represent  the  intersection  of  the  arc  and  the  cone  generated 
by  the  azimuthal  motion  of  the  antenna  (see  Figure  1).  The 
north  and  south  cross  sections  of  the  arc  have  a  maximum 
density  of  1.6  and  0.8  x  10s  cm  '  respectively.  The  angle 
between  geographic  north  and  the  arc  alignment  was  equal  to 
-23.5°.  The  LOS  velocity,  in  spite  of  being  measured,  is  not 
presented  here  because  of  its  poor  statistical  accuracy  and 
its  coarse  spatial  resolution  following  from  the  rapid  radar 
scan. 

Figure  2d  shows  the  slow  AZ  scan  that  starts  at  0230.  The 
LOS  velocity  measured  in  the  F  region  (second  range  gate) 
shows  smaller  values  as  the  counterclockwise  AZ  scan 
progresses.  At  -90°  azimuth  it  changes  from  a  toward 
orientation  to  a  direction  away  from  the  radar.  The  LOS 
velocity  corresponding  to  the  E  region  (first  gate)  shows  a 
similar  characteristic,  with  the  reversal  occurring  earlier  in 
the  scan.  The  statistical  error  of  the  I  OS  velocity  measure¬ 
ment  is  20  m  s  1  for  gate  I  and  about  50  m  s  1  for  gate  2.  In 
spite  of  the  latter  being  a  sizable  error  when  the  magnitude  of 
the  velocity  is  small,  the  velocity  reversal  can  be  visualized 
following  the  trend  of  the  velocity  values.  At  the  end  of  the 
scan  (southwest),  gates  1  and  2  exhibit  LOS  velocities 
oriented  away  from  the  radar  and  with  a  value  of  order  100 
ms1.  The  density  contours  indicated  that  the  arc  was 
oriented  -30°  away  from  the  north. 

The  second  partial  EL  scan,  in  Figure  2c.  measured  the 
complete  cross  section  of  the  arc.  and  also  probed  the 
dawnside  region  outside  the  are  (right  side).  The  region  of 
maximum  density  was  found  displaced  84  km  toward  mag- 
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TABLE  1.  Motion  of  the  S-A  Are  Along  the  Dawn-Dusk  Meridian 


Type  of 
Scan 

Maximum 
Pedersen 
Conductivity, 
x  104  mho/m 

Time. 

UT 

Altitude. 

km 

Dawn-Dusk 

Displacement. 

km 

Time. 

Difference.* 

s 

Distance 

Difference.* 

km 

Velocity.  Dawn 
to  Dusk.+  ms’ 

EL 

1.96 

0225:26 

129.9 

107.6 

AZ 

1.86 

0226:45 

125.3 

152.6 

79 

45.0 

569  (469) 

AZ 

2.43 

0229:50 

120.7 

151. 1 

185 

-1.5 

-8  (-108) 

AZ 

1.95 

0231:37 

131.1 

171.9 

107 

20.8 

194  (94) 

EL 

2.20 

0233:19 

121.2 

170.2 

102 

-1.7 

-161-116) 

AZ 

3.60 

0235:07 

125.2 

175.3 

108 

5.1 

47  (-53) 

AZ 

2.89 

0238:12 

128.5 

187.3 

185 

12.0 

65  (-35) 

AZ 

2.84 

0240:59 

128.5 

188.6 

167 

1.3 

8  (-92) 

•Differences  are  relative  to  the  preceding  scan. 

+  Left-hand  number  Earth  corotating  frame,  right-hand  number  Earth-Sun  frame. 


netic  west  from  its  previous  location  9  min  earlier  (a  dawn  to 
dusk  virtual  velocity  by  the  Nc  contours  of  about  150  m  s  ' 1 ). 
The  LOS  plasma  velocity  within  (and  nominally  transverse 
to)  the  arc  was  about  150  ms1  in  the  F  region  and  300  ms'1 
for  the  E  region,  a  factor  of  3  increase  with  respect  to  the 
velocities  of  Figure  2a.  The  F  region  LOS  velocities  in  the 
dawnside  of  the  Sun-aligned  arc  (right  side  of  the  plot)  are  the 
reverse  of  on  the  duskside,  with  a  value  equal  to  100  m  s  '1. 

For  altitudes  below  150  km  in  the  (production-recombina¬ 
tion)  chemical  equilibrium  region  of  the  Nt.  profile  (exami¬ 
nation  of  the  slope  of,  for  example,  the  0.4-0. 8  x  10 5  cm' ? 
Ne  contours  of  Figures  2 a  and  2e  and  Plate  1)  we  find  the 
duskside  slope  entering  the  arc  shallower  than  the  dawnside 
slope  leaving  the  arc.  This  suggests  a  gradual  hardening  of 
the  precipitating  electron  flux  across  the  arc  from  the  leading 
edge  to  the  trailing  dawnside.  where  the  greatest  depth  of 
penetration  occurs.  There  is  a  relatively  sharp  cutoff  of 
ionizing  electron  flux  on  the  (trailing)  dawnside  of  the  arc. 
We  cannot  claim  here  a  statistically  significant  sample  (three 
out  of  the  three  cases  of  0224,  0233,  0210  on  February  26, 
1987,  and  consistency  for  the  fourth  case  of  Figure  4c  at 
0340,  March  I.  1987). 

The  difference  in  character  (bifurcated  versus  continuous) 
between  the  Ne  contours  of  Figures  2c  and  2g  is  due  to  the 
geometry  of  the  ISR  conical  scan  intersecting  the  slab  of  arc 
ionization.  This  may  be  appreciated  by  visualizing  a  verti¬ 
cally  elevated  slab  or  small  horizontal  cylinder,  either  pierc¬ 
ing  the  (55°  vertex  angle)  vertical  conical  surface  and  passing 
through  near  its  axis,  or  grazing  (and  encompassing  a  portion 
of)  the  conical  surface  to  one  side.  Thus  the  intersection  of  a 
conical  scan  with  Sun-aligned  cylindrical  forms  will  depend 
on  the  distance  of  the  axis  of  the  horizontal  cylinder  from  the 
vertical  axis  of  the  cone,  and  will  look  very  much  like 
Figures  2<7.  If,  2g,  and  2/i.  The  intersection  of  a  conical  scan 
with  a  horizontally  stratified  ionosphere  is  a  set  of  concentric 
circle  Ne  contours.  Figures  2b  and  2c  show  Sun-aligned  arc 
contours  merging  to  the  southwest  with  some  horizontally 
stratified  ionospheric  plasma.  The  core  ionization  in  Figure  2 
(a().8  x  1()'  cm  ’)  is  continuous.  Sun-aligned,  and  quite 
stable  in  density,  size,  and  height  (e.g..  Figures  2 a  and  2c). 
and  separates  F  region  plasma  flowing  away  from  the  Sun  on 
its  dawnside  from  plasma  flowing  much  more  slowly  and 
toward  the  Sun  on  its  duskside  (the  plasma  is  predominantly 
antisunward  within  this  arc). 

Several  characteristics  of  the  measured  LOS  velocity  and 
the  number  density  associated  with  the  arc  of  February  26 


can  also  be  explained  by  changes  in  the  relative  location  of 
the  arc  with  respect  to  the  radar.  The  small  LOS  velocity  of 
Figure  2d  at  -90°  azimuth  is  interpreted  as  a  near  perpen¬ 
dicularity  between  the  flow  and  the  sight  direction.  The 
differences  in  the  peak  density  observed  on  consecutive  AZ 
scans  or  during  two  sectional  cuts  of  the  same  arc  (Figure  2c) 
are  due  to  variations  in  the  distance  between  the  radar  and 
the  S-A  arc.  E  region  velocities  collected  during  the  AZ 
scans  can  be  used  to  reassure  that  the  null  LOS  velocity  is 
due  to  a  near  perpendicularity  of  the  flow  and  the  sight 
direction.  The  velocity  in  the  E  layer,  being  smaller  and 
rotated  clockwise  with  respect  to  the  F  region  flow,  will  give 
a  null  LOS  velocity  at  a  different  look  angle,  as  in  Figure  2d. 

Figure  3  shows  the  density  contours  of  Figures  2d  and  2 f 
superimposed  to  the  Ol  (6300  A)  emissions  (hatched  region) 
recorded  simultaneously  by  the  GL  all-sky  imaging  photom¬ 
eter.  Both  radar  and  optical  data  sets  have  been  transformed 
and  scaled  to  a  common  coordinate  system  that  depicts 
geographic  latitude  and  longitude.  The  red  line  optical  aurora 
was  assumed  to  originate  at  200  km  of  altitude.  The  emission 
produced  by  impact  excitation  shows  a  bright,  discrete,  and 
isolated  feature  located  westward  of  the  radar  station.  The 
number  density  contours  measured  by  the  ISR  mostly  over¬ 
lap  the  hatched  region,  except  for  a  small  part  to  the  south. 
This  region  may  correspond  to  the  oval  F  layer.  Figure  3 
confirms  genuinely  and  unambiguously  that  the  radar,  in 
fact,  diagnosed  an  arc  which  was  Sun-aligntd  and  locally 
produced. 

2.2.  February  26,  1987  Dwells 

Plates  lu  and  16  show  observations  of  the  same  stable  and 
energetic  arc  as  in  Figure  2.  obtained  between  0204  and  0222 
as  the  arc  passed  overhead  the  antenna  parked  in  the 
direction  up  B.  the  magnetic  field.  An  integration  time  of  1 
min  was  selected  in  order  to  provide  temperature  error  bars 
smaller  than  109 f  at  a  200-km  altitude.  The  arc  motion  across 
the  radar  beam  cannot  be  assumed  constant.  However  from 
Figures  2a-2h  (see  Table  I)  we  know  that  the  Sun-aligned 
arc  moved  westward  (toward  dusk)  at  order  of  100  m  s  1  in 
the  Earth  corotating  frame.  As  a  result  of  this  motion,  the 
duskside  of  the  arc  appears  earlier  (lefthand  side)  in  Plates  lu 
and  16.  Color-coded  vertical  profiles  of  T,  and  7,  <  7,  are 
superimposed  on  the  density  contours  in  Plates  lu  and  16. 
respectively.  In  the  center  of  the  arc  the  ion  temperature  is 
500  K  in  the  F  layer  and  about  700  K  at  F  region  heights.  7, 
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Fig.  .3.  Schematic  representation  of  the  auroral  luminosity  vibserved  by  the  GL'ASIP  on  February  2b.  1987.  The 
number  density  measured  by  the  Sondrestrom  (SR  is  superimposed  as  a  thick  line  contour.  The  two  contour  levels 
correspond  to  0.4  and  1.2  x  10 <  cm  '.  Both  6300  A  intensities  and  radar  .V,  have  been  transformed  to  geographic 
coordinates.  They  are  both  colocated  and  aligned  in  the  same  direction  (about  24'  west  of  north). 


sharply  increases  in  the  dawnside  of  the  arc  reaching  1.300  K 
and  exceeding  Te  at  both  E  and  F  layers.  The  electron 
temperature  has  a  different  pattern,  with  T,  enhancements  of 
500  K  colocated  with  regions  of  higher  bottomside  electron 
densities  (more  intense  precipitation).  This  behavior  is  re¬ 
flected  in  Plate  I  h.  where  the  7,  IT,  ratio  is  observed  to 
decrease  from  >  1 .5  to  0.8  near  the  dawn  edge  boundary.  The 
S-A  arc.  defined  by  the  0.8  x  105  contour,  drifted  across  the 
radar  beam  in  10  min.  If  we  take  the  average  westward  speed 
to  be  of  order  100  ms1  in  the  Earth  frame,  w  hich  implies  a 
width  of  about  60  km  for  the  S-A  arc.  we  get  reasonable 
agreement  with  the  0.8  x  10'  contours  of  Figure  2.  The 
maximum  density  was  1.9  x  l()'  cm  1  at  130  km  of  altitude. 
The  width  of  the  arc  and  the  number  density  are  both 
consistent  with  the  view  that  the  arc  properties  are  rather 
stable  over  the  time  period  covered  by  Figure  2  and  Plate  1 
(over  half  an  hour). 

2.3.  March  I,  1 987 

To  demonstrate  that  the  conclusions  to  be  reached  from 
the  February  26  data  do  not  represent  a  unique  circum¬ 
stance.  we  present  data  for  another  Sun-aligned  arc  of 
comparable  intensity,  observed  three  nights  later.  On  this 
night  a  system  formed  by  two  Sun-aligned  arcs  traversed  the 


Sondrestrom  field  of  view  while  the  radar  performed  a 
20-min  sequence  of  AZ-EL-EL  scans.  (At  this  time  the  A.',, 
and  the  solar  flux  indices  were  2~  and  72.8  respectively). 

Figure  4a  shows  the  number  density  and  the  LOS  velocity 
measured  with  the  full  circular  AZ  scan  of  0324  with  the 
antenna  at  a  fixed  elevation  of  45°.  Most  of  the  patchy 
features  to  the  east  and  south  of  the  station  are  radar 
intersections  of  the  polar  cap  F  region.  The  elongated  trace 
located  135  km  westward  of  the  radar  is  the  dear  signature  of 
a  S-A  arc.  We  have  already  noted  that  any  significant 
ionization  enhancement  near  150  km  is  due  to  ongoing 
ionization,  and  will  have  attendant  impact-excited  optical 
emissions.  Thus  we  must  say  something  about  the  oval 
electron  density  enhancements  approximately  80  km  east  of 
the  northern  and  southern  extremes  of  this  elongated  trace. 
These  have  electron  density,  electron  temperature,  and 
electrodynamics  consistent  with  Sun-aligned  arc  signatures 
(piercing  rather  than  grazing  the  conical  radar  scan).  How¬ 
ever.  the  ISR  alone  cannot  ambiguously  identify  that  they 
are  connected  by  continuous  Sun-aligned  enhanced  electron 
density  contours,  and  thus  they  will  not  be  discussed  further 
in  this  paper. 

A  pair  of  EL  scans  (Figures  4  h  and  4c)  cut  this  arc  system, 
respectively  in  its  northwest  and  southwest  regions.  8  and  14 
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Plate  1  ri .  Color-coded  representation  of  the  ion  temperature  as  a  function  of  altitude  and  time  'or  February  26. 
1987,  time  interval  from  0209  to  0221 .  The  i\\.  contours,  measured  simultaneously  superimposed.  The  antenna  was 
pointing  parallel  to  the  local  B  lines  (EL  =  80’.  A 7.  -  141  I.  Each  individual  profile  is  the  result  of  l-min  integration  time. 
Density  contours  start  at  0.2  x  10s  cm  '. 


min  after  it  was  mapped  by  the  AZ  scan.  An  ionization 
feature  to  the  northeast  is  seen  at  0343  on  the  EL  scan  at  51° 
azimuth.  This  (soft  ionizing  flux)  feature  is  not  a  stable 
Sun-aligned  arc.  as  it  does  not  extend  into  the  southeast 
radar  scan.  This  will  be  discussed  later. 

T  he  right  side  of  Figure  4b  (southern  distances)  shows  the 
LOS  velocity  at  160  km  of  altitude  (gate  2)  varying  as 
expected  if  the  velocity  is  constant  and  the  variation  is 
produced  only  by  cnanges  in  the  radar  look  angles.  A 
different  behavior  is  observed  in  the  northern  part  (left  side) 
of  the  plot,  where  the  LOS  velocities  increase  abruptly  and 
then  start  decreasing  at  -100  km  of  distance.  The  region  of 
larger  LOS  velocity  coincides  with  the  dawnward  edge  of  the 
arc . 

3.  Data  Analysis 

T  his  section  presents  the  analysis  of  the  data  correspond¬ 
ing  to  February  26  and  March  1 ,  1987.  Here,  we  describe  the 
method  (hat  was  used  to  calculate  the  velocity  vector,  the 
Hall  and  Pedersen  conductivities,  the  electric  field,  the 
horizontal  and  the  field-aligned  currents,  the  vertical  com¬ 
ponent  of  the  Poynting  flux  and  (he  ion-neutral  differential 
temperature.  T  he  data  acquired  during  the  FT.  and  AZ  scans 


have  been  processed  independently  using  different  algo¬ 
rithms  (  ’ctions  3.1  and  3.2)  which  were  tailored  respec¬ 
tively  for  each  of  these  two  types  of  scans.  The  results  of 
these  independent  calculations  provide  mutually  consistent 
values  of  the  geophysical  parameters,  lending  further  confi¬ 
dence  to  the  analysis  procedures  and  results. 

3.1.  Analysis  of  'Elevation  Scans 

At  F  altitudes  the  plasma  drifts  in  the  ExB  direction.  In 
the  F  layer  the  large  ion-neutral  collision  frequency  prevents 
the  ions  from  moving  perpendicular  to  E.  Brckkc  ct  at.  1 1973. 
I974|  used  the  relation  between  winds,  drifts  and  electric 
fields  to  obtain  the  E  region  neutral  winds.  They  employed 
ion  velocity  vectors  measured  at  two  diffe  ent  altitudes  be 
connected  bv  the  same  field  line. 

A  different  approach  was  exercised  by  Doupnik  el  al. 

1 1 977 )  and  dv  la  Beaujurdiere  el  al.  1 1 977]  in  which  simulta¬ 
neous  LOS  velocities  measured  from  the  E  and  F  layers  and 
the  radar-deduced  wind  velocity  were  used  to  determine  two 
components  of  the  velocity  vector;  the  third  component 
(parallel  to  B)  was  assumed  to  be  null. 

Nevertheless,  the  ion  motion  along  B  is  seldom  zero 
Mechanisms  such  as  adiabatic  expansions,  parallel  electric 
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Plate  l/=.  Similar  to  Plate  In.  hut  the  color-coded  diagram  represents  the  temperature  ratio  i/(  /  l.  The  tune  is  the 

same  a'  in  the  previous  figure. 


fields  or  the  small  component  of  the  neutral  wind  along  B.  all 
drive  ion  motions  along  magnetic  lines.  Span  cr  cl  ul.  [  IV76| 
detected  vertical  neutral  w  inds  in  the  auroral  zone,  less  than 
or  of  the  order  of  80  ms  1 .  Pctchcrxch  cl  ul.  [  1085)  observed 
values  of  15  ms  1  associated  with  auroral  arcs.  Here  we 
measure  the  magnitude  of  the  ionic  motion  in  the  direction 
parallel  to  B  hy  pointing  the  antenna  along  "  on  February  26. 
1087  between  0204  and  0222.  before  the  sequence  of  FI.  and 
AZ  scans  and  during  the  overhead  transit  of  the  arc  \  was 
found  equal  to  25  e  7  m  s  1  In  general  a  nonzero  V  will 
introduce  an  error  equal  to  V  x  tan  (El.)  Consequently  the 
null  l  assumption  fails  when  the  radar  is  directed  near  the 
up  H  position. 

The  velocity  vectors  were  derived  using  the  algorithm 
described  a  Appendix  A.  Phis  method  has  some  of  the 
characteiistics  of  the  one  employed  by  Ihuipnik  cl  ul.  \  1477). 
with  the  main  ditfen_.  ee  that  the  input  data  were  collected  by 
El.  scans  instead  of  dwells.  Appendix  A  also  includes  a 
description  of  the  pulse  conjuration  that  u.ts  used,  and  a 
list  of  the  different  constraints  which  cuitail  the  capability  of 
the  radar  to  measure  the  ionospheric  parameters. 

Concurrently  with  the  radar,  a  Fabry -Perot  interferometer 
recorded  the  intensity  and  the  Doppler  velocity  of  ( »l  infotl 
At  emissions  at  the  zenith  and  lour  cardinal  points  (.1 
Meriwether,  personal  communication.  |dX8i.  I  he  neutral 
wind  pattern  measured  on  I  ebrtiarv  26.  |>)8“  has  K^n 
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tig  -la  Niimtsei  dci.Mlv  and  I  I  IS  velocity  col  responding  so 
March  I.  tug'  the  lime  interval  is  between  and  IH»1  I  he 
antenna  scanned  Ual  m  ,'6oiit  ~  nun.  die  elevation  was  maintained 
lived  at -l'  I  lie  Nun  dnevlion  is  .ndic.ued  ill  the  lovvei  tell  coinci  ol 
die  plol 
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Fig.  4 b.  Elevation  scan  performed  in  the  plane  of  the  magnetic  meridian  between  1.431  and  0336.  Magnetic  north  is 
at  the  left  side  of  the  plot.  Only  arc  A  was  probed  by  the  radar  during  this  scan. 


described  by  Niciejewski  et  al.  [1989],  They  showed  that  the 
direction  of  the  wind  was  southerly  during  the  zenithal 
transit  of  the  arc. 

For  all  the  data  presented  here,  we  have  performed  a 
complete  velocity  vector  analysis,  including  actual  measured 
neutral  winds,  and  a  three-dimensional  magnetic  field  model. 
The  results  of  the  velocity  analysis  for  the  two  EL  scans  of 
Figure  2  are  displayed  in  Figures  5c  and  6c.  Here  we  present 
the  resolved  vectors  and  their  error  bars  (dashed  line). 
Geographic  north  is  at  the  top  of  the  panel  and  the  arrow  in 
the  upper  part  points  toward  the  Sun  location.  The  top 
panels  of  these  figures  show  the  height-integrated  Pedersen 


and  Hall  conductivities  plotted  as  a  function  of  westward 
distance  from  the  radar.  Both  Hall  and  Pedersen  conductiv¬ 
ities  were  calculated  using  the  radar-measured  number  den¬ 
sity.  the  Si  hunk  and  Walker  [1970]  expression  of  the  ion- 
neutral  collision  frequency,  and  the  Banks  and  Kockarts 
I I973J  expression  of  the  electron-neuiral  collision  frequency. 
The  neutral  densities  are  from  the  MS1S-86  model  fHedin. 
1987). 

Figure  5  corresponds  to  the  EL  scan  of  0224.  During  this 
time  the  Pedersen  conductivity  peaked  at  a  horizontal  dis¬ 
tance  of  105  km  away  from  the  radar,  the  Hall  conductivity 
at  95  km.  The  three  vectors  farthest  to  the  west  have  very 
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Fig  4r.  Same  as  Figure  4 b.  but  for  the  elevation  scan  in  a  plane  perpendicular  to  the  magnetic  meridian.  The  cross 
sections  of  both  arcs  are  displayed  at  -  180  and  150  km.  The  polar  cap  F  region  is  on  the  right  side  of  the  plot  at  altitude 
of  300  km 
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large  uncertainties  and  are  given  least  statistical  weighting. 
On  the  contrary.  the  velocities  closer  to  the  arc  have  small 
error  bars  because  of  the  higher  V  These  vectors  show  a 
very  consistent  amisunward  How.  I  lie  ntagnittkle  of  the 
velocity  vanes  from  less  than  ItHI  m  s  1  ;tt  the  center  of  the 
arc  up  to  9(H)  tn  s  1  m  the  limits  ol  the  HI.  scan.  The  region 
of  large  conducliv  Hies  coincides  w  ith  the  small  l  •  1(10  nts  1 1 
velocities.  The  sharp  decline  of  the  conductivities  in  the 
daw  nside  t right  side i  of  the  arc  is  partially  due  tv)  our  limited 
HI .  scan,  which  did  not  piobe  the  /  region  at  distances  close 
to  the  radar. 

figure  o  shows  the  velocity  vectors  from  the  analysis  of 
figure  2 e  ut2'2).  At  this  time  the  region  of  maximum 
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Pedersen  conductivity  was  located  at  175  km  ground  range 
from  the  station  and  the  Hall  conductivity  peak  at  a  horizon¬ 
tal  distance  of  155  km.  Similar  to  figure  5.  the  velocity 
shows  a  steady  antisunward  llow.  Besides,  the  region  of 
v  elocity  gradients  is  extended  in  the  daw  nside  of  the  arc 
reaching  a  value  of  14X4  m  s  1  and  then,  slightly  decreasing 
when  the  Pedersen  conductivity  returns  to  its  background 
level  of  0.2  mhos.  Our  results  agree  with  those  of  Robinson 
<-/  at.  1 1 0X7 1  and  Mcr.dc  el  al.  |I9XX|. 

The  unchanging  characteristic  of  both  Pedersen  and  Hall 
conductivities  during  two  HH  scans  separated  by  X  nun 
argues  in  favor  of  this  S-A  arc  being  very  stable  with  no 
appreciable  variation  of  its  source  and  dynamics.  Based  on 
the  steadiness  of  this  arc.  we  proceeded  to  merge  the  data  of 
f  igures  5  and  h  and  form  the  three  upper  panels  of  figure  7. 
This  new  set  of  v  elocity  v  ectors  and  conductiv  ities  w  as  used 
to  compute  the  height-integrated  and  field-aligned  currents 
following  the  expressions  given  by  llrekke  el  id.  11974]  and 
de  la  Heaiijtiriliere  el  al.  |  1 977 1 . 


J 

-  Iff 

*  -//f 

xBff  -  1/1  x  B 

•  l„il  x  B)  x  B  B 

(1 

./  -v  j 

(2 

where  .1  and  /  are  the  current  perpendicular  and  parallel  to 
B.  1;.  and  2.,/  are  the  height-integrated  Pedersen  and  Hall 
conductivities. 

The  geophysical  parameters  and  their  error  bars  shown  in 
figures  7</-7.g  were  all  calculated  in  a  coordinate  system  for 
which  the  •  .V  axis  is  perpendicular  to  the  tire  alignment  and 
directed  toward  dawn.  The  •  )  axis  is  parallel  to  the  arc  and 
points  sunward  (poleward).  I  he  orientation  of  the  arc  was 
determined  from  the  A/,  scans  displayed  ir  figures  2h-2li.  It 
was  found  that  the  tingle  between  the  arc  alignment  and 
geographic  north  is  equal  to  -24  i  5.  I  he  direction  of  the 
sun  is  7  west  of  north. 

figure  ~d  shows  the  component  of  the  electric  field 
perpendicular  to  the  are  /.  ,  (component  of  the  v  elocity  along 
the  arc).  is  less  than  12  mV  m  1  in  (he  duskvvard  edge  of 
the  arc.  It  reaches  its  minimum  value  of  5.4  mV  in  1  in  the 
region  of  maximum  conductivity  .  Then,  it  gradually  in¬ 
creases  in  the  daw  nside  until  it  presents  a  peak  value  of  72 
mV  m  1  at  94  km  eastward  from  the  center  of  the  arc.  I:\  . 
plotted  in  figure  7c.  has  a  different  behavior.  It  remains 
constant  and  below  6  mV  m  1  except  fora  region  (distant  73 
km  from  the  center)  in  the  extreme  dawnward  end  of  the  arc. 
where  I.  increases  to  2X  mV  ni  1 . 

Both  components  of  the  height-integrated  horizontal  cur¬ 
rents  are  plotted  in  figures  7/  and  7g.  7,  (transverse  com¬ 
ponent)  is  proportional  to  /:, .  it  is  directed  toward  dusk  and 
has  a  maximum  value  of  (l.09h  A  ill  ’.  ./,  .  the  component 
along  the  alignment,  presents  two  different  regions.  Dawn- 
ward  from  the  center  of  the  arc  the  average  current  is  0.05X 
A  ni  1 .  Duskvvard. ./,  is  below  the  uncertainly  level  old. (MIX 
A  m  1 . 

Considering  that  the  K  held  and  the  conductivities  do  not 
vary  along  the  arc  alignment.  /  is  given  by 

a  2.  r  ,12.,!  >'(/-. ,  a/., 

7  - /  -  If" .  i/,  —  Ml 

a  i  a  v  a  i  .i  \ 

Ihe  parallel  current  associated  with  the  arc  of  I  ebruary 
2<).  19X7  was  calculated  following  (4>.  I  lie  numerical  deriv- 
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Fig.  7.  Geophysical  parameters  which  were  derived  using  the 
conductivities  and  velocity  vectors  of  both  El.  scans  of  February  26. 
19X7  (Figures  2 a  and  2c).  la)  Pedersen  conductivity:  l/>)  Hall 
conductivity  :  U  )  velocity  vectors:  (</)  v  component  of  the  K  field  and 
the  statistical  error  inherent  to  the  measurement:  (r)  same  as  Figure 
Id  for  v  component:  If)  r  component  of  the  height-integrated 
horizontal  current  .  Igl  v  component  of  the  height-integrated  hori¬ 
zontal  current:  (It)  field-aligned  current  density.  (/)  ion-neutral 
differential  temperature  computed  from  the  measured  ion  and  neu¬ 
tral  velocity  ;  and  (/I  vertical  component  of  the  Poynting  tins 


atives  were  obtained  by  fitting  a  second-order  polynomial  to 
four  consecutive  values  of  £,  and  In  Figure  7/t.  down¬ 
ward  currents  arc  plotted  as  positive  numbers.  The  second 
and  fourth  terms  of  the  right  side  of  (3)  produce  very  small 
contributions  to j}l  and  can  be  neglected.  In  the  duskside.  the 
parallel  current  is  less  than  I  grA  m  In  the  region  of 
maximum  conductivity.. /[;  is  directed  upward  with  an  aver¬ 
age  value  equal  to  1  /uA  m  :.  Between  156  and  140  km  j; 
decreases  to  a  level  near  zero.  This  decrease  is  due  to  a  large 
I  2  ft  A  m  )  current  originated  by  a  negative  gradient  of  the 
Pedersen  conductivity,  first  term  in  the  right-hand  side  of  (3). 

The  temperature  enhancements  of  the  ion  gas  observed  on 
the  dawnxide  of  the  arc  are  colocated  w  ith  regions  of  large 
plasma  flow.  This  fact  led  us  to  calculate  the  frictional 
heating  due  to  the  ions  being  dragged  through  the  neutral  gas 
by  an  electric  field,  tedder  mid  Banks  ( 1972]  have  indicated 
that  ion  temperature  effects  can  be  described  by  the  ion 
energy  equation.  In  the  steady  state  T,  depends  on  the  vector 
difference  between  ion  and  neutral  velocities  [ Baron  and 
Wand.  19X3].  Hence,  if  ions  transfer  momentum  to  the 
neutrals  sufficiently  rapidly  that  the  neutral  particles  ap¬ 
proach  the  ion  velocity,  only  modest  T ,  enhancements  will 
occur.  Following  the  derivation  of  Baron  and  Wand  |I983] 
for  momentum  transfer,  we  found  that  for  the  number 
density  of  the  S-A  arcs  analyzed  here  the  time  constant  is 
several  hours,  so  for  large  ion  velocities,  large  T,  enhance¬ 
ments  and  heating  rates  should  occur.  We  have  used  the 
radar  ion  velocity  of  Figure  le  and  the  neutral  velocity  from 
the  FPI  to  calculate  the  ion-neutral  differentia!  temperature 
show  n  in  Figure  7/.  Our  calculations  show  that  between  the 
dusk  and  the  center  part  of  the  arc  the  differential  tempera¬ 
ture  is  less  than  100  K.  Substantial  T,  enhancements  occur 
on  the  dawn  edge,  up  to  1600  K  here.  This  is  simply  because 
this  is  where  the  largest  ion  velocities  occur,  while  the 
neutral  particles  do  not  have  time  to  experience  significant 
acceleration. 

Poynting's  theorem  states  that  the  rate  of  energy  flow  per 
unit  area  is  equal  to  E  x  5//.  Here,  we  present  the  calcula¬ 
tion  of  the  vertical  component  of  the  Poynting  flux  (P.) 
which  was  computed  using  the  electric  field  and  the  horizon¬ 
tal  current  obtained  solely  from  radar  measurements.  Fol¬ 
lowing  Ampere's  law  and  the  current  continuity  equation  (2) 
and  assuming  that  (I)  only  the  field-aligned  current  sheets 
contribute  to  fi//.  and  (2)  the  magnitude  of  the  total  dow  n¬ 
ward./  current  is  equal  to  the  magnitude  of  the  upward  jv 
(they  are  considered  to  be  closed  by  horizontal  Pedersen 
currents),  the  P  is  nonzero  within  the  arc  and  given  by 

Jx  -  Ux))  (4) 

where  <./,)  is  the  average  value  of  the  horizontal  current 
measured  at  the  dawn  and  dusk  boundaries  of  the  S-A  arc. 
Figure  7/  displays  the  values  of  P.  calculated  according  to 
this  simple  algebraic  relation.  Negative  vales  refer  to  a 
downward  energy  flux.  The  region  of  more  intense  precipi¬ 
tation  has  a  small  Poynting  flux.  <  I  erg  cm  ■  s  1 .  However. 
I’,  is  maximum,  about  5  erg  cm  :  s  1  further  east  of  the  arc. 
on  the  downside  and  coincident  w  ith  the  region  of  maximum 

,  field.  Consistency  with  these  findings  can  also  be  claimed 
from  less  direct  calculations  (K.  J.  Weber,  private  commu¬ 
nication.  19X9)  w  hich  lead  to  similar  findings  of  a  downw  ard 
Poynting  flux  on  the  dawnxide  of  a  Sun-aligned  arc  diag¬ 
nosed  by  a  rocket  flight  |  We  her  el  al..  [  I989|.  The  time  span 
for  that  segment  of  rocket  data,  however,  was  short  com- 
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pared  to  an  Alfven  bounce  period,  so  temporal  variation 
cannot  be  excluded  front  those  data,  whereas  the  ISR  data 
here  show  steady  conditions. 

The  arc  cross  sections  of  Figure  2  were  also  used  to 
determine  the  motion  of  the  S-A  arc.  The  enhanced  electron 
density  profiles  measured  by  the  radar  were  transformed  to 
Pedersen  conductivity  profiles.  The  maximum  irp  value 
corresponding  to  each  scan  is  tabulated  in  Table  I.  The 
numbers  at  parentheses  in  the  last  column  refer  to  the  arc 
motion  along  a  dawn-dusk  meridian  in  a  frame  fixed  with  the 
Sun-Earth  frame  of  reference  (positive  numbers  represent 
duskward  velocities).  The  arc  velocities  show  some  scatter. 
However,  the  general  trend  is  a  duskward  motion  of  460  m 
s  1  before  0227  UT  followed  bv  a  dawnward  drift  of  60  m 
s 

Figure  8  presents  the  derived  geophysical  parameters  for 
the  EL  scans  of  March  I.  1987.  This  figure  has  the  same 
format  of  Figures  5.  6.  and  7.  However,  when  the  statistical 
uncertainty  of  the  velocity  was  large  or  the  antenna  was 
pointing  close  to  the  up  H  position  the  parameters  were 
deleted  from  the  plot. 

Figure  8  shows  the  derived  parameters  for  the  EL  scan 
that  was  performed  along  the  magnetic  east-west.  Both 
Pedersen  and  Hall  conductivities  present  two  peaks  corre¬ 
sponding  to  both  arcs  of  Figure  4c.  The  background  conduc¬ 
ts  ity  was  0.7  mhos.  Both  arcs  were  located  at  regions  of 
large  velocity  gradients  where  V  •  K  0.  However,  one  of 
them  (arc  B.  see  Figure  4c)  had  a  sunward  flow  and  the  other 
had  ;tn  antisunward  convection.  Figure  8.  as  well  as  Figure 
7.  shows  the  anticorrelation  between  the  magnitude  of  E  and 
the  conductivities.  This  means  that  a  near-zero  E  field  occurs 
when  the  densities  are  the  largest.  The  F  component  of  E  is 
only  -4  mV  m  1  in  the  center  and  dawnside  of  arc  A.  The 
E,  component  varies  from  a  peak  value  of  -90  mV  m  1  in 
the  dawn  edge  of  arc  A  to  *  120  mV  nt  1  in  the  dusk 
boundary  of  B.  J ,  is  proportional  to  with  a  maximum 

value  near  0. 1  A  m  1 .  Jv  reaches  values  up  to  0.075  Am  1 . 
The  field-aligned  current  is  upward  (2  /j. A  m  :)  in  the  region 
of  arc  A:  j;,  is  downward  in  the  duskside  of  arc  B.  reversing 
to  the  highest  upward  value  of  this  plot  (4  fiA  m  ’)  at  the 
center  of  arc  B.  The  ion-neutral  differential  temperature  of 
Figure  8/  and  the  Poynting  flux  of  Figure  8/  show  enhanced 
values  at  the  boundaries  of  arcs  A  and  B. 

7.2.  Analysis  of  Azimuthal  Scans 

Figure  9  presents  the  velocities  which  were  obtained 
following  the  method  outlined  in  Appendix  B.  We  used  data 
from  the  four  partial  AZ  scans  of  Figures  2 h.  2 cl.  If.  and  2 It. 
This  figure  is  in  geographic  coordinates,  where  (he  light 
traces  refer  to  antisunward  convection  and  the  heavy  traces 
indicate  sunward  flow.  The  vectors  corresponding  to  each 
scan  were  shifted  along  the  dawn-dusk  meridian  in  order  to 
compensate  for  any  displacement  of  the  S-A  arc  which 
occurred  between  the  scans.  We  also  averaged  the  velocity 
vectors  when  two  or  more  vectors  were  located  within  ±20 
km  from  each  other. 

The  vector  velocities  of  Figure  9  reproduce  (verify)  many 
of  the  characteristics  that  were  obtained  from  the  analysis  of 
the  EL  scan  data.  Furthermore,  they  present  additional 
features  which  were  not  resolved  by  the  limited  coverage  of 
the  EL  scan.  The  common  features  are  as  follows:  ( I )  There 
is  the  large  antisunward  velocity  in  the  daw  nside  of  the  arc 


o 


SCAN  DISTANCE  (kml 

Fig,  8.  Same  as  Figure  7,  except  (hat  it  corresponds  to  time  03.77 
of  March  l.  1987. 

with  a  peak  value  near  1200  m  s  ~ 1 .  (2)  The  region  of  velocity 
gradients  w  ith  negative  divergence  has  a  width  of  about  175 
km  measured  along  the  dawn-dusk  meridian.  (7)  Antisun¬ 
ward  velocity  decreases  beyond  the  dawn  edge  of  the  arc.  (4) 
There  is  a  high  degree  of  consistency  during  consecutive 
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EAST-WEST  DISTANCE  (Km) 


Fig.  9.  Velocity  vectors  obtained  from  the  four  slow  AZ  scans 
of  February  26,  1987.  A  I  kms*1  vector  is  represented  as  a  segment 
100  km  in  length.  The  inclination  of  the  vector  indicates  the  direction 
of  the  flow. 


antenna  cycles.  The  new  signatures  are  as  follows:  (1)  The 
duskside  of  the  arc  presents  a  strong  sunward  flow,  with  a 
peak  value  of  500  m  s'.  (2)  The  region  of  sheared  flow 
follows  the  arc  alignment  in  the  area  covered  by  the  scan 
(about  700  km).  (3)  The  antisunward  flow  in  the  dawnside 
shows  that  beyond  the  first-order  constancy  of  plasma  flow 
parallel  to  the  arc.  there  is  a  gentle  second-order  variation 
along  the  alignment  of  the  arc.  The  magnitude  of  the  vectors 
decreases  from  1200  ms1  in  the  poleward  part  of  the  scan 
to  900  ms  1  in  the  equatorward  region.  (4)  The  northern  part 
shows  only  a  gradient  but  not  a  reversal.  (5)  Near  the  center 
of  the  plot  the  sunward  flow  shows  a  rotational  reversal 
which  suggests  the  limit  of  the  sunward  motion.  (6)  Further 
duskward  from  the  region  of  sunward  convection,  the  flow  is 
again  antisunward.  There  is  a  sharp  shea-  with  positive 
divergence.  (7)  The  width  of  the  sunward  flowing  region  is  of 
order  100  km. 

3.3.  Analysis  of  Data  Collet  led  With 
the  Antenna  Directed  Up  B 

Derived  parameters  of  the  data  collected  during  the  dwells 
of  February  26.  1987  are  shown  in  Figure  10.  together  with 
the  Poynting  flux  and  Joule  heating  rate  from  the  EL  scans. 
Figures  lOw.  106.  I  Or.  and  I  Or  correspond  respectively  to 
the  Pedersen  conductivity  calculated  from  the  number  den¬ 
sity  obtained  from  the  overhead  transient  of  the  arc.  the 
energy  deposition  from  precipitating  electrons.  Qr.  which 
was  computed  following  the  method  outlined  by  Wichwaret 
al.  [  1975).  the  electron  temperature  measured  at  250  km  and 
the  ion  temperature  from  lit)  km  of  altitude,  f  igure  I0</. 
reproduced  from  Figure  Ij.  displays  the  Poynting  flux  The 
horizontal  scale  was  changed  using  a  constant  arc  motion  •' 
180  ms  '.  Figure  10/  shows  the  Joule  heating  rate,  which 
was  calculated  according  to  the  relation:  .III  - 
i/.(E  +  U  x  B):.  Figure  lOg  presents  the  height-integrated 
ion  cooling  rate  due  to  collisions  with  the  neutrals.  Hanks 
and  Kot  karts  { 1 973 1  showed  that  the  energy  transfer  rate 


from  ions  to  neutrals  can  be  approximated  as  dU,ldt  = 
— /!,7V„ /3,„t 7 ,  -  7„).  where  Nn(N,i  and  T„(T,)  are  the  neutral 
(ion)  density  and  temperature,  and  /3,„  is  the  ion  energy  loss 
rate.  N„  and  T„  were  calculated  using  the  MSIS-86  thermo¬ 
spheric  model  and  an  exospheric  temperature  equal  to  850 
K.  The  upper  and  low<er  traces  correspond  to  the  energy 
transferred  to  N2  particles  (the  major  neutral  constituent  at 
E  region  altitudes)  from  NT  and  07  ions. 

Figures  10«-10r  show  a  prominent  maximum  occurring  at 
0211  UT.  The  coincidence  of  elevated  temperature  and 
enhanced  TLP  indicates  the  intrinsic  relationship  between 
these  parameters.  The  second  peak  in  T may  be  produced 
by  the  decrease  in  the  electron  cooling  rate  due  to  smaller 
values  of  A',,.  Figures  10</-10g  present  a  well-defined  peak  at 
0220.  In  spite  of  the  fact  that  some  of  these  curves  are  not 
completely  independent  (the  circulation  of  some  of  them 
includes  common  input  parameters),  the  similarity  in  their 
shapes  and  location  of  their  maxima  suggests  a  direct  rela¬ 
tion. 

A  more  quantitative  view  of  the  7j.  and  T,  enhancements 
associated  with  the  center  and  dawnside  of  the  February  26 
arc  is  presented  in  Figure  1 1 .  The  four  profiles,  labeled  E  or 
I.  were  obtained  during  the  overhead  transit  of  the  arc.  when 
the  antenna  was  pointing  antiparallel  to  B.  The  dashed  lines 


Fig.  10.  Directly  measured  and  derived  parameters  of  the  Sun- 
aligned  are  of  February  26.  1987:  (u)  height-integrated  Pedersen 
conductivity:  (hi  energy  deposition  from  precipitating  electrons:  (i  t 
7,  at  250  km  altitude.  (i/I  Poynting  flux,  (el  T ,  at  110  km.  (f\  Joule 
heating  rate,  and  let  energy  transferred  from  ions  to  neutrals. 
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Fig.  1 1 .  Electron  and  ion  temperatures  measured  during  the  up 
B  position.  Dashed  lines  correspond  to  0-1 1  UT  and  the  thick  solid 
line  to  measurements  at  0220  UT.  The  thin  trace  represents  the 
neutral  temperature  calculated  using  the  MSIS-86  model. 


are  the  electron  (E)  and  ion  (I)  temperature  profiles  mea¬ 
sured  between  021 1  and  0212  and  w  hen  the  region  of  higher 
density  was  being  probed.  The  thick  continuous  traces  are 
the  7j,  and  T,  profiles  measured  at  0220  and  near  the 
downward  edge  of  the  arc.  The  Tc  curve  corresponding  to 
the  center  of  the  Sun-aligned  arc  shows  an  enhancement  of 
>1000  K  in  the  F  region.  Notice  that  this  profile  increases 
linearly  at  altitudes  above  150  km.  The  T,  curve  measured  in 
the  downside  was  enhanced  by  7(H)  K  in  the  E  and  lower  F 
regions. 

3.4.  Thermal  Analysis  of  the  March  I.  !9H7 
Elevation  Sean 

The  isodensity  contours  on  Figure  4c  are  reproduced  in 
Plate  2.  Superimposed  on  the  /V,  traces  are  two  color-coded 
contours  corresponding  to  the  ion  (blue)  and  electron  (red) 
temperatures.  The  electron  isothermal  contour  of  1000  K 
penetrates  to  lower  altitudes  near  the  regions  of  enhanced 
number  density.  The  ion  temperature  is  enhanced  in  the 
dawnside  of  the  arc  labeled  A.  reaching  a  value  of  2000  K  at 
150  km.  The  duskside  of  arc  B  is  also  a  region  of  high  7, 
values.  T  his  plot  agrees  and  extends  the  observations  show  n 
in  Figures  10  and  II.  T,  and  /,  are  both  enhanced  but  in 
different  regions  which  do  not  necessarily  colocate.  /',  is 
enhanced  at  the  dawn  or  dusk  (or  both)  sides  of  the  Sun- 
aligned  arc. 

Figure  12  has  the  same  format  as  Figure  10.  T  he  curves 
shown  in  (he  seven  panels  were  calculated  based  on  the  IT 
scan  data  of  March  I.  1987.  They  correspond  to  arc  A  of 
Figure  4c.  In  spile  of  the  large  variability  in  the  parameters 
displayed  in  this  figure,  (he  trend  is  the  same:  Regions  of 
enhanced  /',  arc  located  near  the  center  of  the  arc.  and  7, 
enhancements  are  situated  at  the  arc  boundaries. 

4.  Disc  t  SSION 

We  have  extensively  analyzed  two  Sun-aligned  arcs  diag¬ 
nosed  by  the  incoherent  scatter  radar.  Excepting  unusual 
circumstances,  this  requires  an  ISR  observing  mode  tailored 
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to  nominal  temporal  and  geometric  properties  of  Sun-aligned 
ai  - 

The  two  arcs  studied  here  are  relatively  intense  Sun- 
aligned  arcs.  T  his  optimizes  the  statistical  error  bars  on  the 
parameters  derived  from  the  ISR  data.  One  of  the  two  arcs 
has  been  reported  elsewhere  | Xiciejewski  et  at..  I989J.  with 
detailed  analysis  of  its  optical,  ion  production-loss  rates,  and 
ionizing  flux  properties.  Observed  in  the  UV  by  Polar  Bear, 
and  the  visible  from  DMSP  and  ground-based  imagers  with 
time  continuity,  this  is  a  well-documented  example  of  a 
strong  Sun-aligned  arc.  We  show  in  Figure  3  that  the  position 
of  it  s  optical  signature  is  colocated  w  ith  the  ISR  signature.  A 
second  Sun-aligned  arc  is  presented  primarily  to  demon¬ 
strate  that  the  key  findings  here  are  repeatable. 

The  observations  reported  here  apply  directly  to  Sun- 
aligned  arcs  bright  enough  to  be  readily  seen  by  conventional 
all-sky  cameras  or  visible  and  UV  satellite  images.  The 
findings  thus  apply  to  stable  Sun-aligned  arcs  as  reported 
prior  to  the  1980s.  which  are  found  in  the  polar  cap  of  the 
order  of  5'T  of  the  time.  The  physical  processes  that  we 
demonstrate  are  important  in  these  arcs  are.  however, 
candidates  for  applying  to  a  much  broader  class  of  arcs.  The 
degree  to  which  they  represent  the  character  of  the  much 
more  common  weaker  Sun-aligned  arcs  [  Weber  and  Buchan. 
1981:  Carlson  et  at..  1988|.  found  a  third  to  a  half  of  the  time 
in  the  polar  cap.  is  the  subject  of  a  separate  study  now  in 
progress. 

4.1.  ISR  Observing  Mode 

We  have  presented  and  illustrated  an  ISR  mode  which.  we 
maintain,  can  in  general  identify  and  diagnose  Sun-aligned 
arcs  in  the  polar  cap.  Based  on  nominal  properties  of 
Sun-aligned  arcs  found  with  image-intensified  all-sky  pho¬ 
tometers.  the  mode  maps  mesoscaie  (order  l()’  x  It)1  km) 
areas  with  adequate  time  resolution  to  track  the  passage  of 
arcs . 

Arc  identification  is  based  on  a  two-step  rationale.  The 
ISR  signature  of  the  Sun-aligned  arc  is  a  Sun-aligned  sheet  or 
cylinder  of  persistent  enhanced  ionization  below  200  km  (at 
"Fj  region"  or  "photochemical"  altitudes).  For  ionization 
to  persist  much  below  2(H)  km.  it  must  represent  the  presence 
of  ongoing  ionization,  because  its  chemical  lifetime  is  so 
short.  However,  a  region  of  stable  or  persistent  production 
of  Sun-aligned  ionization  must  by  the  same  token  also  be  a 
region  of  stable  or  persistent  Sun-aligned  production  of 
(electron  impact  excited)  optical  emissions.  Thus  when 
looking  at  contours  of  greater  than  10'  electrons  cm  ’  near 
150  km.  with  chemical  lifetime  of  order  of  a  minute,  and  with 
persistence  over  10  min.  vine  cannot  help  bill  also  expect  to 
see  enhanced  airglow  within  the  field-aligned  projection  of 
these  contours.  One  example  of  this  is  shown  in  Figure  3. 


4.2.  Directly  Observed  Barometers 

The  directly  observed  parameters  are  electron  density  , 
electron  anil  ion  temperature,  and  plasma  linc-of-sight  ve¬ 
locities.  T  he  electron  densities  show  persistent  Sun-aligned 
ridges  of  enhancements.  Contours  of  constant  electron  den¬ 
sity  are  Sun-aligned  over  mesoscaie  distances  of  5(MM(XM) 
km  for  altitudes  between  1(H)  and  2(H)  km.  Electron  temper¬ 
atures  are  enhanced  along  magnetic  field  lines  above  these 
enhanced  iidgi  •  of  ionization,  consistent  with  the  expecta- 
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Plalc  2.  Superposition  of  the  contours  corresponding  to  three  different  parameters  measured  by  the  ISR:  electron  density  (black),  Tr  ( red ) 

and  T,  (blue).  The  data  were  taken  on  March  1,  10N7  at  0337:16  UT. 


V  ai  i  xdaki  s  and  C'aki  son  Poi  ah  Arc  Hi  i  c  ikodvn  amic  s  Km  rc;i  i  k  s 


OISTANCE  (km) 

Pig.  12.  Parameters  derived  from  the  left-hand  half  of  the  HI. 
scan  of  0227  UT  on  March  1 .  19X7.  Figures  I2«-I2.e  display  the  same 
geophysical  quantities  as  in  Figure  10.  The  statistical  uncertainties 
of  T,  and  7,  are  25(1  K  and  500  K  respectively. 

tion  of  heating  by  a  (sheet)  flux  of  incident  energetic  elec¬ 
trons. 

On  the  dawnside  of  the  Sun-aligned  arc.  there  is  a  channel 
of  high  ion  temperature,  exceeding  that  of  the  electron  gas. 
This  channel  of  ion  heating  is  where  the  antisunward  plasma 
flow  is  still  near  its  maximum  value.  As  one  moves,  dawn  to 
dusk,  into  the  arc.  the  ion  temperature  falls  to  equilibrium 
with  the  neutral  particle  gas.  and  only  the  electron  gas  has 
significantly  enhanced  temperature.  These  data  confirm, 
refine,  and  extend  the  ion  temperature  enhancement  finding 
and  semiquantitative  interpretation  by  Carlson  cl  al.  |  I9X4|. 
Antisunward  velocity  gradients  arc  both  coaligned  with  the 
Sun-aligned  contours  of  enhanced  (below  200  km  altitude) 
electron  density.  That  is.  where  the  antisunward  plasma  How 
velocity  decreases  in  going  from  dawn  toward  dusk,  ioniza¬ 
tion  is  enhanced:  when  the  dawn-to-dusk  plasma  velocity 
gradient  has  the  opposite  sense  (increases  dawn  to  dusk)  or 
is  constant,  there  is  no  enhanced  ionization  below  2(H)  km. 
These  confirm  a  simple  arc  electrodynamics  interpretation 
(convergent  electric  fields  produce  converging  Pedersen 
currents  whose  continuity  is  maintained  across  the  arc  by 
incident  energetic  electrons). 

4..V  Derived  Parameters 

The  derived  parameters  tire  horizontal  and  Birkeland 
currents.  K  field,  and  Poynting  and  particle  llux.  We  esti¬ 
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mated  the  steady  state  currents  within  the  tire  from  the 
plasma  density  and  velocity  data.  The  Birkeland  currents 
were  calculated  by  two  independent  means.  First  ’he  elec¬ 
tron  density  contours  of  a  cross  section  of  the  arc  (found  to 
match  the  density  contours  several  hundred  kilometers 
upstream  and  dow  nstream  and  on  earlier  and  later  measure¬ 
ments)  were  used  to  estimate  ’he  steady  state  production 
rate  of  the  arc.  and  hence  the  upward  current  above  the  arc 
carried  by  the  estimated  incident  energetic  electron  flux.  In 
the  second  method  the  horizontal  electric  field  gradient 
across  the  arc  was  derived  from  the  plasma  velocity  gradient 
across  the  arc.  This  electric  field,  applied  to  the  ionospheric 
conductivity  (derived  from  observed  electron  density  pro¬ 
files  and  a  standard  model  atmosphere),  led  to  derived 
Pedersen  and  Hall  currents.  Based  on  arc  symmetry  up¬ 
stream  and  downstream,  transarc  gradients  in  calculated 
Pedersen  currents  were  then  used  to  derive  (from  the  hori¬ 
zontal  divergence  of  currents  across  the  arc)  the  variation  of 
Birkeland  currents  across  the  arc.  The  upward  Birkeland 
currents  derived  from  these  two  different  approaches  were 
comparable,  and  both  peaked  near  the  center  of  the  arc  at 
about  1  |U.A  m  :.  Because  downward  Birkeland  currents 
(presumably  carried  by  upgoing  thermal  electrons)  can  be 
derived  only  where  electron  densities  are  large  enough  to 
give  a  useful  1SR  signal,  complete  current  mapping  cannot 
be  done.  However,  this  approach  gave  a  rather  comprehen¬ 
sive  estimate  of  the  ionospheric  current  system  in  and  near 
the  stable  Sun-aligned  arcs  observed,  an  estimate  that  should 
extrapolate  well  to  such  arcs  in  general.  This  is  summarized 
in  Figure  14. 

We  have  performed  three  independent  calculations  of  the 
energy  into  the  Sun-aligned  tire  by  nonparticle  heating. 
While  this  nonparticle  energy  input  is  distributed  across  the 
full  arc.  it  is  strongest  into  the  dawn  edge.  The  three 
calculations  are  motivated  by  recognition  that  energy  exter¬ 
nal  to  the  ionosphere  (ultimately  mechanical  energy  in  the 
solar  wind)  is  carried  down  into  the  ionosphere  along  mag¬ 
netic  field  lines  as  electromagnetic  energy  or  Poynting  flux. 
This  is  dissipated  as  frictional  heating  of  the  ions  dragged 
through  the  neutral  atmosphere,  or  Joule  heating.  Heat, 
initially  deposited  mostly  in  the  ion  gas.  is  rapidly  passed 
onto  the  neutral  gas. 

The  first  calculation  follow  s  the  rationale  of  a  f>K  x  SB 
Poynting  flux  calculation.  We  note  that  the  Sun-aligned  arc 
changes  much  more  slow  ly  along  than  across  its  axis,  and  is 
very  stable  over  times  large  compared  to  a n  Alfven  bounce 
period.  Neglecting  the  spat ial  and  time  derivatives  then,  we 
take  a  cross  section  transverse  to  the  tire,  and  slice  it  into  a 
series  of  adjacent  differential  elements.  We  look  at  the 
electric  field  differential  (dawn  to  dusk  reduction)  and  Ped¬ 
ersen  current  differential  across  each  element  (the  electric 
field  is  in  this  measured  neutral  wind  rest  frame).  I  he  current 
is  then  crossed  against  the  electric  field  strength  to  derive  the 
first  estimate  of  energy,  plotted  in  Figures  7/.  I ()</.  and  12(7. 
and  labeled  Poynting  llux.  Relative  to  an  idealized  measure¬ 
ment  of  Poynting  llux  this  calculation  has  its  shortcomings. 
Yet  we  feel  it  is  useful,  partly  to  underscore  the  real  energy 
source,  and  partly  because  of  the  good  quantitative  agree¬ 
ment  between  the  ergs  cm  :  s  1  found  from  this  calculation 
and  from  (he  following  two. 

The  second  calculation  is  based  on  the  additional  measure¬ 
ment  of  the  ion  gas  temperature,  (liven  the  measured  value 
of  the  plasma  velocity  anil  the  neutral  atmospheric  v  elocity  . 
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we  can  calculate  the  rale  at  which  heat  should  be  going  into 
the  ion  gas.  due  to  the  relative  motion  of  the  ions  through  the 
neutral  gas.  This  leads  to  the  second  estimate  of  energy  flow, 
and  is  verified  by  the  good  agreement  between  the  calculated 
and  the  actually  observed  ion  temperature. 

The  third  calculation  deals  with  the  heat  loss  from  the  ion 
gas.  to  the  neutrals.  It  is  based  on  the  measured  ion  gas 
temperature  and  the  neutral  gas  temperature.  Because  of  the 
cooling  rate  cross  section  dependence  on  ion  composition, 
which  we  do  not  measure,  we  show  two  values  for  the 
extremes  of  100%  NO'  versus  100%  (),  .  Within  this 
relatively  large  uncertainty,  this  third  energy  estimate  agrees 
with  the  first  two  (with  better  agreement  for  the  more  likely 
dominant  ion). 

We  then  applied  tin  approximate  energy  balance  calcula¬ 
tion  to  these  data.  It  was  reasoned  not  only  that  was  it  logical 
to  expect  a  magnetospheric  energy  source  to  drive  the 
high-speed  ionospheric  plasma  flow  through  thermospheric 
frictional-drag  "loading."  but  that  this  energy  dissipation 
should  ultimately  end  up  mostly  in  the  neutral  atmosphere 
below  the  enhanced  ionospheric  densities  in  the  arc.  It  was 
possible  to  estimate  this  latter  quantity  in  two  ways:  ( 1 )  from 
the  steady  state  rate  of  heat  going  into  the  ion  gas;  and  (2) 
from  the  steady  state  rate  of  heat  lost  from  the  ion  gas.  The 
former  is  derived  from  the  conversion  of  ordered  ion  gas 
energy  due  to  coherent  motion  under  the  action  of  the 
driving  electric  field,  to  disordered  energy  (heat)  as  the  ions 
suffer  randomizing  collisions  with  neutral  particles.  This  is 
proportional  to  the  square  of  the  difference  between  the  ion 
and  the  neutral  gas  hulk  velocity  ( Smbhe  and  Chandra. 
1971 1.  The  hitter  is  derived  from  the  observed  difference 
between  the  neutral  atmospheric  and  ion  gas  temperatures. 
Uncertainty  in  the  ion  composition  leads  to  greater  uncer¬ 
tainty  in  this  latter  estimate  than  in  the  former.  However, 
both  of  these  methods  give  a  value  for  the  energy  deposition 
rate  into  the  neutral  atmosphere  of  a  few  ergs  cm  :  s  ' . 
comparable  to  the  Poynting  flux  into  this  region.  These  three 
independent  calculations  thus  lead  to  a  simple  but  reassur¬ 
ingly  self-consistent  interpretation  of  the  data. 

It  is  important  to  note  a  critical  step  in  the  above  argument 
before  proceeding.  We  note  that  can  be  derived  in  two 
complementary  independent  ways,  one  (theoretically)  from 
measured  ion  and  neutral  particle  velocities,  the  other  di¬ 
rectly  from  the  radar-observed  spectra.  Our  calculations  of 
the  differential  temperature  based  on  the  ion  and  neutral 
velocities  show  good  agreement  with  the  observed  T,  en¬ 
hancements  (see  Plate  I  a  and  F  igure  7/).  The  Fabry-Perot 
also  measured  enhanced  neutral  temperatures,  possibly  co¬ 
located  with  the  region  of  larger  ion  temperatures  (R.  Nicie- 
jewski.  personal  communication.  1989).  The  observed  ion 
heating  can  be  explained  solely  on  the  basis  of  a  quantitative 
analysis  of  the  large  observed  plasma  flow  V  relative  to  the 
observed  neutral  atmospheric  rest  frame  l\  Since  all  three  of 
these  parameters  are  observed,  not  modeled,  we  are  thus 
able  to  determine,  not  merely  assume,  that  these  ion  tem¬ 
peratures  tire  explained  solely  on  the  basis  of  this  heating 
mechanism  (within  small  statistical  uncertainty).  It  is  this 
use  of  the  directly  observed  /, .  V.  and  T  that  allows  us  to 
justify  calculating  the  rate  of  heat  flowing  into  the  thermo¬ 
sphere  on  the  basis  of  a  steady  slate  ion  thermal  balance 
argument. 

We  then  point  out  that  the  magnitude  of  this  heating  rate  is 
also  of  significance.  The  Poynting  flux  and  energetic  particle 


flux  deposit  energy  to  the  neutral  atmosphere  near  and  above 
120  km  of  several  ergs  cm  '  s  1  over  the  Sun-aligned  arc  of 
width  of  order  100  km.  This  compares  to  about  0.5  erg  cm  : 
s  1  fora  global  mean  thermospheric  FUV  heating  rate.  Parts 
of  the  arc  experienced  10  times  this  reference  average,  and 
the  average  across  the  arc  was  three  or  four  limes  this 
reference.  Given  that  the  arc  is  of  order  100  km  wide  and 
persists  for  a  few  hours,  it  is  a  significant  local  heat  source. 
It  is  not  a  negligible  heal  source  even  w  hen  averaged  over  its 
motion  across  the  dark  polar  cap  (of  order  10%  to  tens  of 
percent  of  the  global  mean  EUV  and  thus  comparable  to  a 
presently  missing  polar  thermospheric  heating  component 
|e.g. .  Carlson  and  Cron  lcy.  1989|.)  While  the  February  26 
arc  discussed  is  unusually  intense,  the  polar  cap  under 
northward  IMF  conditions  typically  has  several  stable  Sun- 
aligned  arcs  of  a  tenth  or  more  of  this  visual  (particle 
precipitation)  intensity  over  a  2000  km  or  less  polar  cap 
width,  w  ith  temporal  persistence  of  order  an  hour  or  more. 
Estimates  based  on  measured  particle  precipitation  energy 
flux  alone  under  these  conditions  may  underestimate  the 
thermospheric  heating  rate  due  to  stable  Sun-aligned  arcs  by 
a  factor  of  5  or  more. 

Note  that  for  the  March  I  Sun-aligned  arc  the  particle 
precipitation  energy  (Figure  12)  is  about  a  sixth  that  for  the 
February  26  case  (Figure  10)  although  both  arcs  have  com¬ 
parable  incident  Poynting  flux. 

Here  w  e  have  stressed  taking  advantage  of  the  stability  of 
the  arc.  in  time  and  in  space  along  its  axial  (Earth-Sun) 
dimension,  to  optimally  diagnose  its  character.  Ultimately, 
temporal  and  spatial  variations  must  occur  and  become  large 
enough  to  be  of  interest  during  one  series  of  arc  observa¬ 
tions.  During  the  40  min  observation  of  the  February  26. 
1987  arc.  it  drifted  from  overhead  (0212  UT)  to  100  km  west 
(0225)  to  180  km  west  (0254).  A  composite  image  of  the  arc 
was  constructed  from  the  many  conical  and  planar  cross 
sections  of  the  arc  to  estimate  departures  from  idealized 
constancy  along  its  axial  length.  This  analysis  indicates  a 
velocity  gradient  along  the  arc  axis,  with  penetration  of 
sunward  drifting  plasma  only  part  way  into  the  polar  cap. 
The  two-dimensional  cross  section  (Figure  9)  suggests  flow 
from  where  the  are  connects  to  the  aurora!  oval,  partial  entry 
into  the  polar  cap.  and  (hen  clockw  ise  rotational  flow  back  to 
nntisunward.  exiting  again  from  the  polar  cap.  Farther  into 
the  polar  cap  there  is  slowing  of  antisunw  ard  plasma  flow  , 
but  not  a  reversal  to  sunward.  Phis  is  sketched  schematically 
in  Figure  15.  (This  is  very  reminiscent  of  the  interpretation  of 
some  !)F  data  by  Hoffman  ct  at.  (1985)  under  similar 
conditions.) 

It  is  clear  from  Figure  15  that  the  net  velocity  difference 
from  the  dawnside  to  the  duskside  of  the  arc  is  greater  in  the 
velocity  reversal  region  than  the  merely  velocity  reduction 
region.  Thus  it  is  apparent  from  examination  of  this  figure 
that  the  Pedersen  current  divergence,  and  balancing  Birke- 
land  currents,  lead  to  extra  energetic  particle  input  within 
the  tongue  of  sunward  moving  plasma,  relative  to  that  in  the 
region  of  simply  a  velocity  gradient  deeper  in  the  polar  cap. 
Furthermore,  while  Hall  currents  along  the  dawn  edge  of  an 
are  of  constant  properties  along  its  axial  length  simply 
continue  into  the  polar  cap.  those  in  the  region  of  partial 
penetration  of  sunward  flow  penetrate  only  partially  into  the 
polar  cap  and  then  must  rotate  westward  near  theii  region  of 
maximum  poleward  penetration.  These  effects  can  combine 
to  produce  the  "bright  spot"  or  extra  particle  precipitation 
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Fig.  13.  Schematic  representation  of  the  height-integrated  den¬ 
sity  (upper  side)  and  the  plasma  flow  pattern  (lower  surface).  This 
figure  is  based  on  the  HI.  and  AZ  scans  of  February  26.  (9X7. 


often  seen  near  where  the  Sun-aligned  arc  connects  to  the 
nightside  auroral  oval. 

Finally,  we  condense  the  essential  character  we  have 
found  as  characteristic  of  the  electrodynamic,  thermal,  and 
energetic  properties  of  these  Sun-aligned  arcs  into  the  illus¬ 
tration  of  Figure  14.  The  character  of  these  two  sample  arcs 
(February  26  near  midnight  MLT  and  March  1  post-midnight 
MLT.  1987).  understood  within  the  context  of  Figure  14, 
extends  and  extrapolates  to  at  least  the  broad  class  of  strong 
Sun-aligned  arcs.  We  are  now  combining  ISR  and  ASIP  data 
for  a  large  number  of  stable  Sun-aligned  arcs,  to  establish  the 
extent  to  which  these  findings  characterize  typical  polar  cap 
conditions  for  northward  IMF  conditions  | Carlson  ct  til.. 
1988). 

The  qualitative  features  of  this  "cartoon  "  must  apply  to 
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t  ig.  14  Cross-sectional  view  of  the  Sun-aliened  arc  of  February 
2f>.  J9X7.  F.lectrical  parameters  of  the  are  and  enerj»\  going  into  the 
ionosphere  are  also  indicated.  This  representation  is  in  a  semiquan- 
titative  fashion. 


Sun-aligned  arcs  in  general.  Within  the  center  of  the  arc. 
where  At,,  is  enhanced  below  200  km.  T,  must  be  enhanced 
by  the  incoming  flux  of  impact  ionizing  electrons.  These 
must  also  carry  an  incoming  particle  energy  flux,  and  an 
outgoing  Birkeland  current,  to  one  side,  while  there  must  be 
a  return  (incoming)  current  on  the  other  side.  The  Pedersen 
currents  must  decrease  from  dawn  to  dusk,  as  the  antisun- 
ward  plasma  drifts  must  decrease  from  dawn  to  dusk.  The 
plasma  drifts  may  simply  decrease,  or  reverse  as  in  this 
illustration.  The  ion  temperature  will  be  enhanced  where  the 
difference  between  the  plasma  and  neutral  atmospheric 
velocities  is  sufficiently  great.  This  depends  on  the  thermo¬ 
spheric  w  inds  across  the  polar  cap.  and  the  extent  to  w  hich 
plasma  drifts  may  have  strong  sunward  flow  on  the  duskside 
of  the  arc.  7,  will  generally  be  enhanced  on  the  dawnside. 
and  occasionally  on  the  duskside.  Likewise,  the  Poynting 
flux  typically  peaks  on  the  dawnside,  but  peaks  on  the 
duskside  if  there  are  strong  sunward  flows  on  the  duskside. 
The  magnitude  of  the  Poynting  flux  is  many  ergs  cm  :  s  1 
for  relatively  intense  Sun-aligned  arcs,  exceeding  the  parti¬ 
cle  energy,  and  smaller  by  a  factor  to  be  statistically  defined 
for  more  common  weaker  arcs. 


5.  Conclusions 

This  study  has  led  to  the  following. 

1.  The  first  direct  measurement  over  full  mesoscale 
horizontal  areas  (of  order  ID1  x  10  "  km)  of  persistent  plasma 
velocity  gradient  across  long-lived  ( sr  1 0 ^  seconds)  stable 
Sun-aligned  arcs  hit'-  been  made.  The  sign  of  these  gradients 
confirms  simple  arc  electrodynamics.  This  is  true  for  single 
and  double  Sun-aligned  arcs  observed. 

2.  New  experimental  diagnostic  techniques  have  been 
developed,  including  illustration  of  an  ISR  observing  mode 
capable  of  diagnosing  polar  electrodynamic  features  with 
temporal  stability  of  at  least  several  minutes  and  spatial 
coherence  of  at  least  100  to  a  few  hundred  kilometers: 
extension  of  a  potential  fitting  procedure  to  map  polar 
plasma  flow  from  ISR  data:  and  illustration  of  a  plasma 
rotational  feature  near  where  a  Sun-aligned  arc  merges  with 
the  midnight  auroral  oval. 

3.  The  cross  sectional  area  map  of  the  plasma  density 
and  thermal  character  across  stable  Sun-aligned  tires  has 
been  measured.  The  electron  gas  temperature  is  enhanced 
over  the  arc.  supporting  the  expectation  that  the  arc  is 
produced  by  a  sheet  of  incoming  energetic  electrons. 

4.  It  has  been  determined  that  along  the  high-velocity 
dawnside  of  the  incoming  energetic  electron  sheet  (produc¬ 
ing  the  arc)  there  is  a  channel  of  enhanced  ion  temperature 
within  which  the  ion  temperature  exceeds  that  of  the  elec¬ 
tron  gas.  This  Sun-aligned  channel  within  which  heat  flows 
from  the  ion  gas  to  the  electron  gas  coincides  with  the  high 
antisunward  velocity  channel  along  the  dawn  edge  of  the 
Sun-aligned  region  of  enhanced  electron  density.  In  addi¬ 
tion.  if  there  is  also  return  sunward  flow  at  sufficiently  high 
speed,  ion  heating  is  also  found  on  the  duskside  of  the  arc. 
The  ion  gas  heating  can  be  explained  solely  on  the  basis  of 
ion  frictional  drag  of  high  plasma  velocity  in  the  thermo¬ 
sphere  rest  frame. 

5.  Persistent  strong  Joule  heating  driven  by  Poynting  flux 
down  into  a  Sun-aligned  arc  has  been  calculated  and  found  to 
be  at  a  rate  well  exceeding  the  particle  energy  flux.  This  rate 
is  estimated  by  three  independent  calculations  based  on  a 
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horizontal  current  differential,  a  velocity  differential,  and  a 
temperature  differential.  Furthermore,  the  magnitude  of  this 
newly  measured  heat  flux  is  several  ergs  cm  :  s  1  and  as 
such  is  of  significance  to  the  polar  thermospheric  energy 
budget,  and  contributed  to  solution  of  the  missing  polar 
thermospheric  heat  source. 

6.  Present  knowledge  of  the  basic  character  of  stable 
Sun-aligned  arcs  in  the  polar  cap  ionosphere  has  been 
extended  to  a  relatively  complete  yet  simple  description, 
consistent  with  theoretical  interpretation. 

Appendix  A 

We  show  in  this  appendix  that  the  LOS  ion  velocity 
measured  by  an  ISR  during  elevation  scans  can  be  used  to 
derive  the  full  velocity  vector  (V).  The  only  conditions  for 
the  determination  of  V  are  that  the  LOS  velocity  from  both 
E  and  F  layers  needs  to  be  measured  almost  simultaneously, 
and  the  probed  volumes  need  to  be  connected  by  the  same 
magnetic  line. 

In  the  F  region  the  ion  motion  is  controlled  by 

V  =  E  x  B/fl:  (All 

and  for  the  E  region 

V  =  II  -  1-  A.',U.  x  B IB  +  k,EIB  +  k: E  x  B /B2  (A2) 

where  E.  V.  and  L  are  the  electric  field,  the  ion  velocity 
and  the  neutral  wind  in  a  plane  perpendicular  to  B.  With 

"A 

k  1  -  — ; - ;  (A3) 

"f 4  «r 

=  — ;  (  A4) 

p-  +  n- 

is  the  ion-neutral  collision  frequency  and  if,  the  ion 
gyrofrequency.  The  line  of  sight  velocity  ( V , ,s )  depends  on 
the  velocity  vector  according  to 

V  Eos  —  ^  •  A  (A3) 

where  A  is  the  unit  vector  in  the  antenna  pointing  direction 
expressed  in  magnetic  coordinates. 

The  transmitted  pulse  pattern  consisted  of  two  pulses,  a 
narrow  pulse.  30/as  (4.5  km)  long,  and  a  long  pulse.  320  /as 
(48  km)  in  length.  Data  sampling  time  was  set  equal  to  8  /as 
and  a  32-lag  autocorrelator  was  used  to  provide  the  ac 
function  of  the  incoherent  scatter  signal.  This  setup  allows 
the  spectra  to  accept  contributions  from  a  range  of  82  km. 
The  volume  probed  by  the  single  long  pulse  during  one- 
integration  period  is  equal  to  the  product  of  82  km.  the  arc 
segment  scanned  during  the  integration  time  and  the  radar 
beam  width.  In  the  E  region,  the  number  density,  the  ion 
velocity  and  the  neutral  wind  vary  with  altitude  very  rapidly: 
in  this  sense  the  long  pulse  measurements  are  not  able  to 
resolve  narrow  layers  and  only  provide  a  weighted  average 
over  the  entire  E  region.  This  weighting  factor  depends  on 
the  local  number  density  the  inverse  of  the  range  squared 
and  a  triangular  function  to  compensate  for  the  fact  that  all 
the  lag  products  have  different  spatial  resolution  which 
linearly  decrease  for  the  longer  delays.  F.xpressions  ( A I )  and 
(A2)  were  weighted  following  the  procedure  described  be¬ 
fore.  then  multiplied  by  the  A  vectors  to  form  dot  products 


for  the  E  and  E  layers  \Bnkkt'  <-/  al..  1973)  and  made  equal 
to  the  measured  LOS  velocities. 

Radar  observations  are  subject  to  several  types  of  con¬ 
straints  w  hich  need  to  be  considered  in  any  calculation  of  the 
velocity  vectors.  One  of  them  is  related  to  the  drift  of  the 
S-A  arcs.  I  hese  arcs  move  along  the  dawn-dusk  meridian, 
and  the  Earth  rotational  motion  will  modulate  the  measured 
velocities  resulting  in  the  following  distortions.  (I)  The 
Doppler  velocity  is  reduced  by  the  term  V,.-A.  where  V,.  is 
the  corotation  velocity  .  At  Sondrestrom  V,  is  equal  to  180  m 
s  '.  (2)  There  is  also  a  small  spatial  displacement  between 
contiguous  integration  periods,  equal  to  (Vj  ^  V’.ia. )  x  J, 
w  here  is  the  velocity  of  the  arc  with  respect  to  a  fixed 
frame,  and  T  is  the  time  interval  between  E  and  E  region 
measurements  (about  30  s).  Although  this  effect  is  small  (6 
km),  it  could  generate  significant  errors  when  the  E  region 
measurements  are  mapped  to  the  F  region.  The  large  range 
resolution  (82  km)  hinders  also  the  ability  to  resolve  the 
vector.  In  particular,  low-elevation  measurements  can  in¬ 
clude  a  large  contribution  of  the  horizontal  variations  of  the 
ion  flow  (i.e..  plasma  boundaries). 

in  a  strict  sense  the  E  to  F  region  mapping  algorithm  can 
be  employed  only  for  data  that  have  been  collected  in  a  plane 
coincident  with  the  B  field;  for  Sondrestrom.  this  corre¬ 
sponds  to  -39°  azimuth.  During  our  experiments  the  EL 
scans  were  performed  at  -129°  and  51  azimuth:  conse¬ 
quently.  the  E  region  projections  of  the  E  region  measure¬ 
ments  are  located  slightly  north  and  w  est  of  the  measured  E 
layer  values.  However,  for  the  data  presented  here,  the  arcs 
were  aligned  -24°  and  -25°  with  respect  to  geographic 
north,  and  the  equipotential  contours  (aligned  with  the  arc) 
deviate  14°  or  15’  from  the  magnetic  meridian.  The  E  region 
intersections  were  between  1  or  2  km  away  from  equipoten¬ 
tial  lines  that  crossed  the  E  region  measurements. 

In  order  to  gain  some  confidence  in  the  analysis  of  data 
from  EL  scans,  we  performed  a  zero  order  calculation  of  the 
velocity  vector  using  the  LOS  velocities  measured  at  1 1 Oand 
160  km  of  altitude  and  105  km  of  distance  (see  arrow  in 
Figure  2e).  For  simplicity  U  is  considered  to  be  zero  and  B 
to  be  pointing  downward.  After  subtracting  the  corotation 
term  the  radar  LOS  velocities  are  Vlos(E)  =  173  m  s  '. 
Vjos<E)  -  -101  ms  1 .  Equations  (Al)  and  (A5)  indicate 
that  the  velocity  component  along  the  scan  is  -168  ms1. 
Then,  using  (2).  k ,  -  0.292  and  k:  =  0.298:  the  other 
component  (perpendicular  to  the  scan)  is  -968  ms1.  After 
a  rotation  to  geographic  coordinates,  we  obtain  V  =  739 
m  s  1  and  Vvoll(h  =  647  ms  ’.  Instead,  when  we  use  all  the 
terms  of  (A2).  the  1GRF-85  model  of  B.  a  second  Vj()S 
measurement  in  the  E  region  and  the  neutral  wind  from  the 
Fabry-Perot.  the  result  is  \  c.ls,  =  647.  V'siulth  =  900  ms1. 
This  is  a  velocity  that  points  closer  to  the  antisunward 
direction.  All  the  derived  velocities  (V)  presented  in  Figures 
4-8  were  obtained  using  the  more  precise  and  complete 
algorithm  that  considered  all  the  terms  of  (AIMA5).  The 
zero  approximation  shown  above  can  be  employed  when  a 
computer  is  not  available. 

Appendix  B 

The  LOS  velocity  measured  during  radar  azimuthal  scans 
can  also  be  used  to  derive  the  ionospheric  electric  field  E.  on 
the  velocity  vector,  using  the  property  of  E  of  being  irrota- 
tional.  The  electric  field  is  derivable  from  a  scalar  potential 
following 
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E  =  -VV  (Bl) 

where  V  is  the  electrostatic  potential  that  varies  as  a  function 
of  latitude  and  longitude.  Now.  if  the  scan  is  performed  at  a 
low  elevation,  then  V  can  be  represented  by  V’(r.  it).  where  r 
is  range  and  a  the  azimuth. 

An  incoherent  backscatter  radar  measures  the  radial  com¬ 
ponent  of  the  ion  velocity  at  several  distances  along  the  line 
of  sight.  At  high  latitudes  this  radial  drift  is  closely  related  to 
the  azimuthal  component  of  E.  Consequently.  V  could  be 
determined  by  integrating  the  azimuthal  component  of  (Bl) 
and  using  a  know  initial  value  of  the  potential.  VI  r.  </,,).  At 
the  Millstone  Hill  radar,  the  boundary  potential  can  be 
assumed  to  be  null  at  the  southernmost  part  of  the  scan  \Holi 
a  al.,  1 984 J .  However,  no  similar  statement  can  be  used  to 
reduce  data  from  Sondrestrom. 

Another  method  that  has  been  used  to  resolve  the  magni¬ 
tude  and  the  direction  of  the  ion  flow  is  the  beam  swinging 
technique.  Here,  the  LOS  velocity  from  several  adjacent 
azimuthal  sectors  are  least  squared  fitted  to  extract  the  three 
components  of  the  velocity  [Ihifilors  and  Bclinkc.  1 974] . 
This  method  has  been  employed  at  Arecibo.  where  the  more 
uniform  ionosphere  of  this  mid-latitude  station  grants  confi¬ 
dence  for  the  application  of  this  type  of  analysis. 

Although  Hr.  «„)  could  be  calculated  from  the  velocity 
vectors  of  section  3.1.  we  instead  computed  the  boundary 
potential  using  the  data  from  the  same  azimuthal  scan.  We 
emploved  the  LOS  velocity  of  three  adjacent  integration 
periods  and  the  condition  of  negligible  ion  motion  along  B 
(less  than  or  of  order  25  m  s  1  as  measured)  to  obtain 
space-averaged  velocity  vectors.  Then,  the  initial  boundary 
potential  was  calculated  for  every  integration  period  and 
served  as  the  initial  condition  for  the  more  precise  value  of  V 
which  was  obtained  by  integrating  (Bl).  The  boundary 
potential  as  calculated  in  this  way  smooths  rapid  spatial 
variations  of  the  velocity  but  will  be  less  affected  by  statis¬ 
tical  errors.  This  scheme  was  adopted  to  avoid  cumulative 
errors  produced  by  random  fluctuations  in  Vj()S  which 
sometimes  are  of  the  same  order  of  magnitude  as  the 
measured  velocity. 
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